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SECTION  I 
INTRODUCTION 

1 . BACKGROUND 

Two  basic  criteria  must  be  satisfied  to  provide  accurate  design  or  analysis 
of  structures.  First,  the  applied  loading  conditions  on  the  structure  must  be 
defined  and,  second,  structural  behavior  models  must  be  available  which  adequate- 
ly simulate  the  structure's  response  under  the  applied  loading  conditions.  The 
calculations  involved  can  range  anywhere  from  simple  hand  solutions  to  complex 
computer  codes.  The  required  complexity  of  the  loading  function  will  depend 
upon  the  degree  of  refinement  attempted  in  analyzing  the  structure. 

For  an  aboveground  structure  subjected  to  airblast  loading,  the  external 
boundary  conditions  arc  composed  of  (a)  the  airblast  induced  pressures  acting 
on  the  exposed  face;  of  the  structure  and  (b)  the  foundation  reactions  resulting 
from  these  pressure  loadings  and  the  motions  of  the  structure. 

These  airblast  loadings  have  generally  been  determined  by  one  of  the  follow- 
ing methods. 

a.  Extension  of  experimental  data  derived  from  high  explosive  tests  and 
shock  tube  and  wind  tunnel  experiments. 

b.  Theoretical  calculations  based  on  ideal  gases. 

c.  Large  hydrodynamic  codes  employing  theory  and  experimental  data  as  the 
basis  of  their  behavior  models. 

Since  most  experiments  have  been  conducted  at  incident  overpressures  less  than 
50  psi  and  only  a few  tests  above  100  psi,  the  use  of  empirical  methods  without 
verification  becomes  questionable  at  higher  pressures;  extrapolation  can  be 
misleading  and  ideal  gas  behavior  may  no  longer  apply. 

2.  OBJECTIVE 


The  primary  objective  of  this  investigation  was  to  determine  the  adequacy  of 
current  techniques  for  predicting  airblast  loadings  on  exposed  aboveground 
structures  subjected  to  high  incident  overpressures.  To  satisfy  this  objective, 
experimental  data  were  collected  to  evaluate  and  refine  current  simplified  tech- 
niques and  complex  mathematical  models.  Because  the  initiaV’diffraction  loading 
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phase  was  considered  the  area  of  greatest  uncertainty,  this  investigation  was 
directed  at  the  early  time  loading  of  simple  structural  shapes.  Two  basic 
shapes,  a rectangular  block  and  a cylinder,  were  subjected  to  nominal  100  and 
300  psi  incident  overpressure  levels.  Pressure  measurements  were  recorded  on 
the  front,  roof,  and  rear  section  of  the  rectangle  and  at  specified  intervals 
around  the  upper  front  and  rear  of  the  cylinder.  Additional  measurements  were 
taken  on  a simulated  ground  plane  in  front  of  both  structures.  The  effect  of  a 
horizontal  roof  extension  on  both  structures  was  also  investigated. 

The  experimental  data  were  used  to  provide  both  a description  of  the 
phenomenology  associated  with  high  pressure  shock  diffraction  and  an  evaluation 
of  current  simplified  prediction  techniques  specified  by  The  Air  Force  Manual 
for  the  Design  and  Analysis  of  Hardened  Structures.  An  evaluation  of  a more 
complex  analytic  technique  for  predicting  airblast  diffraction  is  reported  in 
reference  1 and  employs  the  data  collected  under  this  effort. 

3.  APPROACH 

The  experimental  phase  consisted  of  testing  basic  two  dimensional  structural 
shapes  in  a 13-inch  diameter  high-pressure  shock  tube.  The  time  duration  of 
interest  and  height  and  depth  of  the  test  structures  were  such  that  the  effects 
of  the  third  dimension  are  not  inherent  in  the  test  presented.  The  shock  tube 
used  in  this  investigation  was  a combustible  gas  driven  tube  located  at  the 
Civil  Engi  'ering  Research  Facility,  Kirtland  AFB,  New  Mexico.  Nominal  100  and 
300  psi  st  shock  pulses  were  used  as  incident  loading  conditions  on  the  test 
structures  shown  in  figure  1.  Representative  waveforms  of  each  loading  are 
shown  in  figure  2.  Figures  l.(a)  and  l.(b)  show  the  basic  rectangular  and 
cylindrical  structures  and  also  indicate  the  central  gage  locations  for  mea- 
suring pressure.  The  supports  for  both  structures  were  designed  such  that  each 
structure  could  be  rotated  between  tests  allowing  measurements  at  all  locations 
around  the  structures.  Figures  1 . (c)  and  l.(d)  show  the  wedges  used  to  slice 
the  incident  shock  wave  and  provide  a surface  for  measuring  reflected  shocks 
propagating  from  the  faces  of  the  structures.  The  orientation  view  of  t.iese 
figures  is  such  that  the  shock  propagation  direction  is  down  and  parallel  to 
the  plane  of  the  figure.  Two  gage  locations  were  chosen  on  the  left  faces  of 
the  wedges.  These  gages  were  on  a line  parallel  to  the  direction  of  shoe 
propagation  and  intersecting  the  stagnation  point  on  the  structure  faces. 
Because  the  gages  lie  along  a plane  of  synmetry  with  respect  to  the  shock  front 
and  structure  face,  the  test  setup  may  be  viewed  as  measurements  along  a ground 
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(a)  Rectangular  Structure 


(c)  Rectangle  and  Wedge  Extension 


(d)  Cylinder  and  Wedge  Extension 


Figure  1.  Test  Structures  (Continued) 
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(e)  Rectangle  with  Roof  Extension 


if)  Cylinder  with  Roof  Extension 


Figure  1.  Test  Structures  (Continued) 
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plane  which  intersects  the  midpoints  of  the  structure  face.  Figures  l.(e)  and 
l.(f)  show  the  roof  extensions  attached  to  the  rectangle  and  cylinder.  These 
were  tested  to  determine  whether  the  pressure  loading  on  the  front  and  top 
surfaces  would  be  independent  of  the  downstream  extensions.  In  addition  to 
the  on-structure  pressure  measurements,  measurements  of  incident  pressure  were 
made  during  each  test  for  control  purposes.  Additional  details  concerning  the 
test  structures  and  gage  locations  are  provided  in  figures  3 and  4. 

After  conversion  to  engineering  units,  the  experimental  data  were  evaluated 
to  produce  details  associated  with  the  shock  diffraction  process--detai Is  con- 
sisting primarily  of  reflected  pressures,  relief  times  to  stagnation  pressure, 
and  vortex  action  on  the  top  and  rear  of  the  structures.  Once  the  basic 
phenomenology  was  described,  the  data  were  compared  to  predictions  based  on 
techniques  specified  in  the  available  literature.  These  comparisons  were  then 
used  to  evaluate  the  validity  of  the  prediction  techniques. 
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Figure  3.  Rectangular  Test  Structure  Details 
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SECTION  II 
TEST  DESCRIPTION 


1.  SHOCK  TUBE 

A general  description  of  the  13-in.  high-pressure  shock  tube  is  shown  in 
figure  5.  Also  depicted  in  this  figure  are  the  locations  of  the  pressure  gages 
used  to  monitor  the  incident  pressure  during  each  test  and  the  velocity  count 
and  oscilloscope  start  triggers.  The  combustible  gas  for  all  tests  consisted 
of  one  part  oxygen  to  two  parts  hydrogen.  The  15-ft.  long  combustion  chamber 
is  sufficient  to  produce  a constant  or  step  function  shockwave  for  approximately 
400  microseconds  at  the  test  station  located  65  feet  away. 

2.  TEST  PREPARATION 

After  sealing  both  ends  of  the  shock  tube  expansion  chamber  with  thin 
aluminum  sheeting,  the  hydrogen  and  oxygen  gases  were  fed  into  the  chamber  and 
allowed  to  mix  for  a minimum  time  of  1 hour.  Prior  to  this  time  the  test 
structures  had  been  instrumented  and  positioned  in  the  test  section.  The  test 
section  was  then  sealed  by  closing  the  access  ports.  The  instrumentation 
recording  systems  were  given  a final  check  and  ambient  measurements  of  tempera- 
ture and  pressure  were  made.  Firing  of  the  explosive  gas  mixture  was  accom- 
plished by  spark  plugs  equally  spaced  along  the  top  and  bottom  lengths  of  the 
combustion  chamber.  The  oscilloscopes  used  to  record  the  transducer  signals 
were  automatically  activated  by  the  trigger  gages  located  upstream  from  the 
test  structures.  These  procedures  were  repeated  for  each  test. 


15  FT  -<— ■ 
COMBUSTION 
CHAMBER 


65FT 

EXPANSION  CHAMBER 


24IN- — 7INH- 

[ VE>  OCITY  COUNT  | J-  SCOPE  START 
j :/»'  SIDE  ON  PRESSURE  GAUGE 


8'/.  i II  IN  ACCESS  PORT  MODEL  SUPPORT 
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Figure  5.  Thirteen-inch  High-Pressure  Shock  Tube 
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3.  TEST  STRUCTURES 

The  test  structures  used  in  this  investigation  are  shown  in  figure  1.  The 
3x3x12  in.  rectangular  structure  was  instrumented  with  pressure  transducers 
located  on  the  upper  half  of  one  face  as  shown  in  figure  3.  It  was  necessary 
to  space  these  gages  along  the  third  dimension  (width)  in  order  to  reduce  the 
vertical  spacing  between  gages.  For  the  short  time  durations  of  interest  and 
the  dimensions  of  the  shock  tube  and  test  structures,  the  variations  in  pressure 
at  each  gage  location  will  be  independent  of  the  third  dimension.  Further, 
symmetry  considerations  require  the  measurement  of  pressure  only  on  one-half  of 
the  vertical  face  of  the  test  structure.  Measurements  were  conducted  on  the 
roof  and  rear  faces  by  rotating  the  test  structure  in  the  manner  suggested  by 
figure  6.  However,  pressure  measurements  were  made  only  at  two  upper  locations 
on  the  rear  face  of  the  model  since  the  proximity  of  the  lower  locations 
to  the  ..model  support  plates  would  have  yielded  questionable  data  from  reflec- 
tions and  f,'ow  disturbance. 


R.  FRONT  FI!!  I.  IMF-FRONT  H I IF  C.  ROOF-IEII  H R IF  0.  NERI 


Figure  6.  Rectangular  Structure  Rotation  Sequence 

The  3-in.  diameter  by  12-in.  length  cylindrical  structure  was  instrumented 
with  three  gages  at  10  degree  intervals  as  shown  in  figure  4.  Again,  the  third 
dimension  was  used  to  allow  closer  angular  spacing  of  the  gages.  The  cylindri- 
cal structure  was  designed  such  that  the  central  4 inch  section  could  be 
rotated,  thus  providing  complete  pressure  measurement  capability  around  the 
structure.  Naturally,  numerous  tests  were  necessary  to  provide  complete  data 
acquisition  around  both  structures. 
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Symmetry  considerations  allow  the  top  half  of  the  cylindrical  structure  to 
be  viewed  as  an  arch  lying  on  a planar  surface  when  boundary  layer  effects  are 
neglected.  For  this  reason,  the  cylindrical  structure  will  be  referred  to  as 
an  arch  throughout  the  remainder  of  this  report. 

All  test  structures  were  supported  across  the  center  of  the  shock  tube  by 
a backup  plate  which  was  bolted  to  either  side  of  the  shock  tube  wall.  This 
plate  also  served  to  conduct  the  gage  cables  from  the  test  structure  to  the  out- 
side of  the  shock  tube.  The  support  spacers  shown  in  figures  3 and  4 were 
replaced  with  the  roof  extensions  when  these  tests  were  conducted. 

4.  INSTRUMENTATION 

Model  ST4  and  ST4A  pressure  transducers  manufactured  by  Susquehanna  Instru- 
ments were  used  to  measure  incident,  stagnation,  and  on-structure  pressure 
loadings.  These  gages  differ  only  in  that  the  ST4  gage  is  designed  for  a 
threaded  receptacle  whereas  the  ST4A  is  not  (see  figure  7).  The  manufacturer's 
specifications  for  these  gages  are  presented  in  appendix  A.  These  gages  are 
miniature  piezoelectric  bar  type  pressure  transducers  which  have  been  found  to 
be  suitable  gages  for  high  pressure  shock  tube  experimentation.  The  ST4  and 
ST4A  gages  have  a high  frequency  response  with  low  acceleration  sensitivity  and 
low  cross-axis  sensitivity. 

At  the  higher  300  psi  test  level  the  gages  had  to  be  protected  against  the 
temperatures  associated  with  the  high  reflected  pressures  and  gas  combustion 
products.  This  protection  was  provided  by  a layer  of  Scotch  electrical  tape 
covered  with  a thin  layer  of  silicone  lubricant  as  an  ablative  material. 

Initial  data  recording  of  these  gages  was  accomplished  with  Kistler  charge 
amplifiers,  Models  503  and  568,  and  Tektronix  555  oscilloscopes  using  1A6 
preamplifiers  and  recorded  on  Polaroid  film. 

Calibration  of  the  pressure  transducers  was  accomplished  with  the  use  of  a 
specially  modified  hydraulic  dead  weight  tester.  This  device  was  capable  of 
providing  step  pulses  of  fixed  magnitude.  Each  gage  was  calibrated  at  various 
pressure  levels  and  sensitivity  settings.  The  resulting  calibration  curves 
were  used  to  convert  the  Polaroid  recordings  of  oscilloscope  traces  to  final 
engineering  units.  Periodic  recalibration  was  performed  during  the  testing 
phase  to  detect  any  changes  in  gage  response  characteristics.  If  and  when 
detected,  the  new  calibration  curve  was  used. 
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5.  TEST  SEQUENCE 

A complete  list  of  all  tests  conducted  during  this  investigation  is  provided 
in  appendix  B.  Included  in  this  list  are  the  ambient  pressures  and  temperatures 
recorded  prior  to  each  test.  In  general,  the  test  sequence  consisted  of 
(1)  preliminary  tests  with  no  structure  present  to  determine  the  correct  amounts 
and  mixture  for  the  explosive  gas  which  would  yield  100  and  300  psi  incider L 
pressures  at  the  test  section  and  (2)  tests  conducted  with  a fully  instrumented 
structure  in  the  test  section.  With  the  use  of  a stagnation  probe  mounted  in 
the  center  of  the  test  section,  the  first  group  of  tests  also  provided  informa- 
tion on  reflected  and  stagnation  pressures  on  the  face  of  the  probe  at  varying 
incident  pressure  levels.  The  stagnation  probe  was  in  the  shape  of  a frustum 
of  a right  cone.  The  tip  of  the  frustum  had  a diameter  only  slightly  larger 
than  the  pressure  transducer  mounted  there  (see  figure  7). 
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SECTION  III 
DATA  REDUCTION 

Thorough  interpretation  of  test  data  often  requi res  an  understanding  of  how 
the  raw  experimental  data  were  treated  during  conversion  to  final  computer 
produced  plots.  This  section  presents  a general  description  of  the  process 
used  to  convert  the  raw  data  to  engineering  units. 

1.  OPTICAL  - TO  - DIGITAL  CONVERSION 

All  gage  response  data,  including  calibrations,  were  initially  recorded  as 
oscilloscope  traces  on  Polaroid  film.  These  traces  were  then  reduced  from  op- 
tical form  to  digital  data  with  a Universal  Telereader  in  the  following  manner. 
The  oscilloscope  traces  were  electronically  converted  to  digital  records  as 
the  operator  of  this  device  manually  followed  each  trace  with  moveable  cross- 
hairs. The  output  from  this  device  in  counts  goes  to  an  X and  Y accumulator 
which  then  transfers  the  data  to  a summary  card  punch.  A calibration  for  each 
trace  was  obtained  by  recording  the  X and  Y coordinates  of  diagonal  points 
on  an  indicated  4x4  -cm  grid  on  the  Polaroid.  During  all  optical  reduction, 
a 10  power  magnification  was  used  resulting  in  a resolution  of  approximately 
400  counts/cm  where  the  data  traces  were  confined  to  a minimum  5 x 10  -cm  grid 
area.  The  resulting  punched  card  output  was  then  stored  on  magnetic  tape  in 
card  image  format. 

2.  DIGITAL  - TO  - ENGINEERINGS  UNITS  CONVERSION 

Before  the  test  data  could  be  converted  to  standard  engineering  units  of 
pressure  and  time,  it  was  necessary  to  generate  calibration  values  for  each 
pressure  transducer.  This  was  accomplished  by  converting  the  digital  analog 
of  the  calibration  oscilloscope  traces  to  units  of  time  and  deflection  based 
on  the  oscilloscope  grid  scale  and  the  prescribed  time  scale  tor  each  calibra- 
tion. A typical  time  duration  for  each  calibration  pulse  was  four  seconds. 

Each  pressure  time  curve  was  then  numerically  integrated.  The  value  of  this 
integral  just  prior  to  pressure  relief  was  then  divided  by  the  time  duration 
to  provide  an  average  deflection  or  amplitude.  The  deflection  values  thus 
generated  and  the  known  calibration  pressure  levels  were  then  used  to  provide 
pressure  versus  deflection  calibration  curves  for  each  transducer.  All  data 
conversion,  integration,  and  plotting  were  performed  by  a computer. 
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After  the  calibration  curves  were  generated,  it  was  then  possible  to  con- 
vert the  test  data  from  digital  format  to  units  of  pressure  and  time.  The  com- 
puter program  developed  for  this  data  reduction  first  converted  the  digital 
data  to  deflection  versus  time  coordinates  and  then  used  the  linear  segmented 
calibration  curves  to  convert  to  units  of  pressutr  versus  time. 

The  label  appearing  on  each  individual  data  plot  is  used  to  identify  the 
plot  according  to  the  designation  system  described  by  figure  8.  The  exception 
to  this  system  is  the  identifying  label  used  for  the  stagnation  probe  tests. 

The  following  system  was  used  for  these  tests: 

1 - STAG.  PROBE  - INC  * 296 
or 

1 - STAG.  PROBE  - FF  * 297 

where  the  "1"  represents  the  test  number,  INC  represents  an  incident  pressure 
measurement,  FF  represents  a measurement  on  the  probe  tip,  and  the  last  char- 
acters represent  the  measurement  number. 

3.  COMPOSITE  DATA  CONSTRUCTION 

When  two  or  more  measurements  were  made  at  similar  locations  and  under 
similar  conditions,  composite  data  records  were  constructed  by  averaging  the 
similar  measurements.  Similar  measurements  were  averaged  only  when  there  were 
no  major  discrepancies  in  waveform  or  magnitude  between  measurements.  The 
composite  record  identification  system  used  for  these  traces  is  the  same  as 
for  the  individual  measurements  with  the  exception  that  all  measurement 
numbers  used  to  generate  the  composite  are  listed. 
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SECTION  IV 
EXPERIMENTAL  DATA 

1.  INTRODUCTION 

The  objective  of  this  section  is  to  provide  characteristic  descriptions  of 
the  test  data  in  a condensed  format.  The  stagnation  probe  data  are  presented 
primarily  as  reflection  factors  and  drag  coefficients  as  functions  of  incident 
overpressure  level.  Test  data  from  the  simple  structural  shapes  are  shown  in 
the  form  of  pressure  variations  as  a function  of  time  and  position.  Local 
reflection  factors  and  drag  coefficients  are  also  presented.  The  data  from  the 
extended  model  tests  are  provided  as  overlays  on  the  pressure-time  histories 
recorded  at  similar  locations  on  the  basic  structural  shapes  without  extensions. 
And, finally,  the  simulated  ground  plane  test  data  are  presented  as  reflected 
wave  attenuation  curves. 


' 


More  specific  information  in  the  form  of  individual  pressure-time  histories 
are  presented  in  appendix  C for  all  measurements  conducted  during  this  effort. 
Appendix  C also  contains  composite  pressure-time  histories  constructed  from 
similar  measurements.  These  composites  were  used  in  the  following  data  pre- 
sentation when  possible. 

2.  STAGNATION  PROBE  TESTS 

As  stated  earlier,  preliminary  tests  were  conducted  using  varying  amounts 
and  ratios  of  oxygen  and  hydrogen.  The  objective  of  these  tests  was  to  identify 
the  correct  gas  mixtures  which  would  yield  the  desired  100  and  300  psi  incident 
pressure  levels.  Table  1 provides  a tabulation  of  incident  pressures  and  re- 
sulting reflected  and  stagnation  pressures  measured  during  these  tests.  In- 
cluded in  this  tabulation  are  the  shock  front  velocities  which  were  measured 
with  trigger  gages  located  upstream  from  the  shock  tube  test  section.  A com- 
parison of  the  measured  incident  pressures  and  incident  pressures  derived 
theoretically  from  the  measured  shock  front  velocities  is  also  provided  in 
table  1.  This  comparison'  demonstrafes  reasonable  agreement  between  the  two 
quantities  and  indicates  that  within  this  range  of  pressures,  free  stream  shock 
characteristics  do  not  vary  significantly  from  ideal  behavior.  In  fact  the 
primary  source  of  variation  may  easily  be  attributed  to  experimental  errors  in 
measuring  pressure  and  shock  velocity. 
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Table  1 

STAGNATION  PROBE  MEASUREMENTS 


INCIDENT  REFLECTED 
PRESSURE  PRESSURE 

(psi)  (psi ) 


quantity  Not  Measured 


STAGNATION 

PRESSURE 

(psi) 


SHOCK  FRONT 
VELOCITY 
(fps) 


100 

510 

275 

104 

570 

310 

104 

585 

280 

104 

650 

375 

106 

530 

290 

107 

675 

350 

110 

650 

350 

113 

650 

350 

135 

745 

390 

153 

950 

550 

158 

920 

460 

161 

980 

465 

224 

1590 

925 

300 

2140 

1250 

330 

2350 

1400 

CALCULATED 
INCIDENT 
PRESSURE 
(psi ) 
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The  data  collected  during  these  tests  can  be  used  to  determine  reflection 
factors  and  drag  coefficients  for  a range  of  incident  pressure  levels.  The 
values  would  be  applicable  to  frontal  loads  on  a disc  shaped  plane  oriented 
normal  to  the  direction  of  shock  propagation.  The  reflection  factors  may  be 
calculated  simply  as  follows: 


RF 


where 


RF  = Reflection  Factor 


= Reflected  Pressure 


P^g  = Incident  Overpressure 


The  drag  coefficients  require  a knowledge  of  the  density  (p)  and  flow  velocity 
(u)  behind  the  shock  front.  Neither  of  these  quantities  were  measured  in  the 
experiments.  However,  p and  u may  be  calculated  from  the  measured  incident 
pressures  and  ambient  shock  tube  conditions  as  follows: 

Co 

n + W7 


5P 


u = 


so 


7P, 


P 


7 + 6VPo 
7 + Pso7Po 

_ 


where 

C = Ambient  Speed  of  Sound 
o 

p = Ambient  Density  of  Unshocked  Gas 
o 


The  above  relationships  are  for  an  ideal  gas  with  constant  specific  heats.  How- 
ever, within  the  range  of  pressures  investigated  there  is  negligible  difference 
in  flow  velocity  when  air  is  considered  an  ideal  gas  or  a nonideal  gas  (refer- 
ence 2)  and  the  variation  in  density  of  the  shocked  gas  amounts  to  less  than 
10  percent  at  300  psi  under  ideal  and  nonideal  considerations. 
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With  values  of  p and  u determined  as  above,  the  dynamic  pressure  may  be 
calculated  as  qo  - 1/2  pu*.  The  difference  between  measured  quasi-steady 
pressure  and  incident  overpressure  is  commonly  defined  as  drag  pressure,  For 
t e stagnation  probe  tests,  this  guasi-steady  pressure  is  the  same  as  the 
stagnation  pressure.  The  experimental  drag  coefficients  were  then  determined 
the  ratio  between  drag  pressure  and  dynamic  pressure. 

Reflection  factors  and  drag  coefficients  determined  from  the  preceding 
relationships  are  presented  in  table  2.  Figure  9 shows  reflection  factors  as 
a function  of  incident  overpressure  level  and  figure  10  displays  similar  infor- 
mation for  drag  coefficients.  Both  sets  of  data  are  weighted  heavily  by  the 
larger  number  of  lower  pressure  level  measurements.  Figure  9 shows  the  expect- 
ed increase  in  reflection  factors  with  increasing  overpressure  level  while 
igure  10  shows  a smaller  Increase  in  drag  coefficients.  A disturbing  observa- 

on  is  the  relatively  larger  scatter  in  data  shown  in  figure  ,0  when  compared 
to  the  reflection  factor  data  of  figure  9. 
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Table  2 

STAGNATION  PROBE  REFLECTION  FACTORS  AND  DRAG  COEFFICIENTS 

I INCIDENT  REFLECTION  DRAG 

I PRESSURE  FACTOR  COEFFICIENT 

(psi) (RF) (Cd) 


94 

4.79 

1.12 

95 

4.89 

1.16 

100 

5.1 

1.32 

104 

5.48 

i.46  ; 

104 

5.63 

1.25 

104 

6.25 

1.93 

106 

5.00 

1.27 

107 

6.31 

1.66 

110 

5.91 

1.57 

113 

5.75 

1.49 

135 

5.52 

1.25 

153 

6.21 

1.64 

158 

5.82 

1.20 

161 

6.09 

1.17 

224 

7.10 

1.76 

300 

7.13 

1.66 

330 

7.12 

1.66 
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INCIDENT  OVERPRESSURE  (psT ) 


Figure  10.  Stagnation  Probe  Tests  Drag  Coefficients 


3.  STRUCTURAL  MODEL  TESTS 

The  experimental  data  gathered  during  the  structural  model  tests  are  pre- 
sented in  three  separate  subsections.  These  subsections  consist  of  the  basic 
structural  shapes,  the  extended  roof  models,  and  the  models  with  wedge  exten- 
sions (ground  planes).  Each  of  these  subsections  is  further  divided  into  the 
data  associated  with  the  arch  and  data  associated  with  the  rectangle. 

a.  Basic  Structural  Shapes 
(1 ) Simple  Arch 

General  descriptions  of  pressure  loading  around  the  arch  structure 
at  nominal  100  and  300  psi  incident  pressure  levels  are  presented  in  figures  11 
and  12.  These  figures  show  the  pressure  distributions  at  12.5  psec  intervals 
between  0 and  400  psec  for  the  100  psi  level  and  between  0 and  300  psec  for  the 
300  psi  level.  Inspection  of  these  figures  provides  a description  of  the  basic 
phenomena  associated  with  shock  diffraction.  First,  envelopment  of  the 
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Figure  11.  100  Psi  Overpressure  Load  Distribute 

Arch  Structure  (Continued) 


INCH  (KIKES) 


Figure  12.  300  Psi  Overpressure  Load  Distribution  - Arch  Structure 
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structure  by  the  incident  shock  wave  is  seen  to  result  in  subsequently  higher 
pressures  on  the  upstream  face  from  reflection  behavior.  As  time  progresses, 
the  reflected  pressure  levels  are  eventually  reduced  to  a quasi-steady  level 
by  rarefaction  waves.  The  resulting  pressures  on  the  windward  face  are  then 
composed  of  the  combined  action  of  overpressure  and  a positive  dynamic  pressure 
effect  (drag  loading).  Both  the  reflected  pressure  levels  and  the  later  quasi- 
steady pressure  levels  are  seen  to  decrease  at  higher  angles  of  incidence. 

The  rear  face  locations  show  that  the  initial  pressure  levels  are 
generally  lower  than  the  incident  overpres'  jre.  This  decrease  results  from  the 
expansion  of  the  incident  pressure  wave  into  the  rear  region  of  the  structure. 
After  the  initial  passage,  the  pressure  levels  continue  to  be  less  than  the 
incident  overpressure.  However,  the  reduction  is  now  caused  by  the  action  of 
vortices  which  follow  the  initial  wave  and  are  generated  as  the  shocked  gas  flow 
deflects  around  the  crown  of  the  structure.  These  vortices  generate  extremely 
complex,  dynamic  loading  behavior  during  the  first  half  of  the  diffraction  pro- 
cess, but  eventually  form  somewhat  stable  hign  and  low  pressure  regions  on  the 
rear  face.  The  loading  throughout  this  period  of  time  is  composed  of  the  shock 
wave  overpressure  and  again  a dynamic  pressure  effect  However,  in  this  situa- 
tion the  dynamic  pressure  results  in  a negative  drag  e net  due  to  the  direction 
of  the  flow  field  which  is  additionally  disturbed  by  the  localized  vortices. 

The  "valleys"  displayed  on  the  curves  suggest  the  location  and  number  of  vor- 
tices acting  at  any  one  time.  In  the  later  stages  of  the  diffraction  process 
the  number  and  magnitude  of  these  vortices  tend  to  stabilize,  but  not  disappear. 
A comparison  between  the  100  psi  and  300  psi  data  indicates  that  the  lower  pres- 
sure level  is  subject  to  a larger  number  of  localized  vortices  than  the  higher 
pressure  level. 

Tables  3 and  4 provide  a list  of  initial  peak  pressures,  steady 
state  pressures,  and  time  lapses  between  initial  peak  load  and  the  steady  state 
load.  Also  included  in  these  tables  are  the  measured  incident  pressures,  mea- 
sured shock  front  velocities,  and  incident  pressures  calculated  from  the  shock 
front  velocities.  Again,  the  latter  quantities  are  provided  as  a measurement 
check.  Since  each  test  on  the  arch  structure  included  three  different  on-struc- 
ture measurements  (different  locations),  the  measured  shock  front  velocity  and 
calculated  incident  pressure  are  listed  only  once  for  each  group  of  measurements 
common  to  any  one  test. 
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Table  3 

ARCH  STRUCTURE  MEASUREMENTS-lOO  PSI  OVERPRESSURE 


INCIDENT 

PRESSURE 

(PSI) 

109 

109 

109 

109 

109 

109 

m 

109 

109 

no 

— no- 


INITIAL 

PEAK 

PRESSURE 

(PSI) 

525 

607 

500 

550 

555 

630 

467 

560 

438 

440 

m 

3o4"~ 


STATE 

PRESSURE 

(PSI)., 

300 

323 

320 

305 

310 

341 

295 

313 

285 

260 


TIME 

LAPSE 

(vsecl 

320 

320 

320 

320 

320 

320 

310 

300 

330 

300 

30(T 

“ITT 


FRONT 

VELOCITY 

(FPS) 

INCIDENT 

PRESSURE 

(PSI) 

3217 

3224 

99 

99 

3213 

99 

3196 

99 

50 

no 

310 

199 

290 

i r\  r\  I 

60 

— m 

240 

175 

290 

3238 

1 00 

70 

105 

110 

no 

112 

T3! 

131 

163 

153 

100 

90 

80 

90 

270 

245 

260 

280 

* 

3206 

98 

” " r\r\ 

80 

108 

no 

no 

1 05 
130 
125 

70 

60 

75 

240  ~ 

220 

175 

3221 

99 

90 

108 

no 

no 

no 

80 

67 

88 

105 

40 

40 

60 

70 

T40 

180 

185 

160 

3ZcO 

99 

100 

105 

no 

87 

80 

80 

90 

285 

305 

110 

no 

45 

50 

260 

i n r\ 

rar- 

no 

52 

70 

145 

3199 

98 

125 

108 

34 

50 

2l5 

3197 

93 

130 

no 

35 

90 

K 235~ 

135 

108 

25 

40 

215 

140 

no 

50 

50 

255 

150 

108“ 

30 

30 

245 

160 

108 

37 

60 

235 

3194 

97 

170 

108 

50 

60 

215 

3208 

98 

To! 

108 

103 

20 

20 

20 

75 

75 

55 

1 55 
155 
140 

T80- 

10! 

35 

— 35"“ 

Tor- 

* Quantity  Not  Measured 
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On  the  upstream  face  of  the  simple  arch  structure  (locations  less 
than  90  degrees),  the  initial  peak  pressure  values  in  tables  3 and  4 may  be 
interpreted  as  reflected  pressure  levels.  At  values  greater  than  90  degrees, 
the  same  column  of  values  represent  the  pressure  levels  of  the  initial  wave 
passage  which  has  been  reduced  below  the  upstream  incident  overpressure  level 
because  of  the  expansion  of  the  incident  shock  as  it  deflects  around  the  rear 
face.  The  pressure-time  waveforms  became  increasingly  complex  on  the  downstream 
face  of  the  test  structure  and  quite  often  failed  to  stabilize  into  what  might 
be  identified  as  a true  steady  state  pressure  level.  For  this  reason  some  of 
the  indicated  steady  state  pressures  listed  in  tables  3 and  4 at  angles  of 
incidence  greater  than  90  degrees  represent  average  pressure  levels  over  a late 
time  interval. 


The  dynamic  nature  of  the  structural  loading  requires  that  the 
time  lapse  between  the  initial  loading  and  the  development  of  steady  state 
loading  becomes  an  important  parameter  in  the  description  of  the  diffraction 
process.  The  time  lapse  values  specified  in  tables  3 and  4 represent  one  of 
two  things  depending  on  the  location  of  the  measurement.  In  general,  for 
locations  on  the  front  face  of  the  structure,  the  time  lapse  values  indicate 
the  time  required  to  reduce  the  initial  peak  reflected  pressure  to  the  steady 
state  level  by  the  action  of  rarefaction  waves.  At  higher  angles  of  incidence 
the  time  lapse  values  represent  an  approximate  time  required  for  the  lower 
initial  peak  pressure  to  rise  to  the  average  steady  state  level.  The  time 
lapses  represent  the  combined  effect  of  the  expansion  process  of  the  upstream 
incident  overpressure  into  the  rear  region  of  the  structure  and  th_'  stabiliza- 
tion of  vortex  behavior  on  the  rear  face. 


Initial  peak  pressure  ratios  and  drag  coefficients  were  determined 
from  the  data  in  tables  3 and  4 and  are  presented  in  tables  5 and  6.  The 
initial  peak  pressure  ratios  are  defined  as  the  peak  initial  pressure  measure- 
ments divided  by  the  incident  overpressure  common  to  that  measurement.  The  drag 
coefficients  were  determined  in  the  same  manner  as  described  earlier  for  the 
stagnation  probe  tests.  The  pressure  ratios  for  values  of  theta  less  than  90 
degrees  represent  reflection  factors.  Pressure  ratios  and  drag  coefficients  are 
shown  graphically  in  figures  13  and  14  as  functions  of  theta.  Evident  in  figure 
13  are  the  two  regions  containing  the  extreme  angles  for  Mach  reflection.  The 
angle  of  incidence  at  which  the  two  possible  angles  of  reflection  based  on 
theory  are  equal  is  termed  the  extreme  angle.  Mach  reflection  occurs  when  the 
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Table  6 

ARCH  STRUCTURE  PRESSURE  RATIOS  AND  DRAG  COE FF I C I ENTS- 300  PSI  OVERPRESSURE 


INCIDENT 

Eirag 

THETA 

PRESSURE 

INITIAL  PEAK 

COEFFICIENT 

(DEGREES) 

(PSI) 

PRESSURE  RATIO 

(C„> 

0 

319 

6.58 

1.51 

303 

5.74 

1.29 

10 

319 

5.72 

1.59 

319 

6.11 

1.30 

20 

284* 

7.57* 

1 .42* 

303 

5.78 

1.20 

30 

284*  “ 

5.88* 

1.10* 

303 

5.35 

0.90 

40 

284*'“ 

4.91* 

0.89* 

315 

5.00 

0.80 

50 

315 

4.70 

0.72 

60 

315 

2.63 

0.55 

70 

304 

1.64 

0.16 

311 

2.22 

0.07 

330 

1.09 

-0.03 

80 

305 

1.13 

0.07 

293 

1.21 

0.11 

~90~ 

330 

0.98 

0.02 

293 

0.99 

-0.17 

100 

330 

0.76 

-0.12 

110 

300 

0.75 

-0.09 

120 

300 

0.42 

-0.31 

300 

0.42 

-0.31 

125 

333 

0.38 

-0.11 

135 

333 

0.30 

-0.05 

140 

283 

0.35 

-0.16 

145 

333 

0.11 

-0.20 

150 

283 

0.15 

-0.34 

340 

0.12 

-0.21 

170 

293 

0.51 

-0.08 

340 

0.26 

-0.17 

180 

340 

or 

-0.17 

293 

0.34 

-0.30 

* Questionable  Data 
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Figure  13.  Initial  Peak  Pressure  Ratios  - Arch  Structure 
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angle  of  incidence  is  slightly  greater  than  the  extreme  angle.  The  100  psi  data 
indicate  an  angular  region  located  near  the  40  degree  location  while  the  300 
psi  level  shows  a slightly  greater  angular  position.  This  higher  position  is 
to  be  expected  based  on  theoretical  values  of  the  extreme  angle  which  tend  to 
increase  at  higher  shock  strengths.  However,  this  theoretical  increase  is  small 
compared  to  the  change  indicated  in  figure  13.  The  larger  change  may  be  attrib- 
uted to  a combination  of  two  factors.  First  and  most  important,  the  theoretical 
values  are  determined  for  a planar  surface  intersecting  the  shock  path  and 
second,  any  variations  arising  from  experimental  error  may  tend  to  exaggerate 
the  difference. 


The  drag  coefficients  shown  in  figure  14  indicate  little  pressure 
level  dependency  at  angles  of  incidence  less  than  60  degrees.  At  larger  angles 
of  incidence  the  effects  of  pressure  level  seem  more  pronounced  with  slightly 
higher  values  of  drag  coefficients  indicated  for  the  300  psi  level.  The  data 
at  larger  angles  of  incidence  also  indicate  different  variations  of  drag  coef- 
ficients with  theta  between  the  two  pressure  levels.  The  difference  results 
from  the  number  and  locations  of  the  vortices  acting  on  the  rear  surface  for 
the  100  psi  tests  and  the  300  psi  tests. 


The  time  lapse  values  shown  in  table  3 for  the  100  psi  arch  tests  show 
that  the  period  of  reflected  pressure  relief  remains  essentially  constant  for 
angles  of  incidence  between  zero  and  80  degrees.  At  80  degrees  and  above,  the 
time  required  to  achieve  steady  state  conditions  deci eases  with  a rapid  decrease 
at  90  degrees.  This  location  of  decrease  is  to  be  expected  because  this  is  a 
region  of  small  reflections  and  since  this  region  is  closest  to  lower  free 
stream  pressures,  it  is  consequently  closer  to  the  origin  of  rarefaction  waves. 
On  the  rear  of  the  arch  structure  the  time  lapse  values  indicate  stabilization 
times  for  the  vortices  acting  on  the  rear  surfaces.  Again,  over  a certain  re- 
gion the  time  lapse  values  are  relatively  constant;  the  region  being  between 
90  and  160  degrees.  Below  this  region  a significant  decrease  is  noted.  For 
the  300  psi  data  shown  in  table  4 the  regions  of  essentially  constant  time  lapse 
are  found  between  0 and  60  degrees  and  between  100  and  170  degrees.  The  re- 
maining areas  are  then  affected  by  more  rapid  changes.  The  obvious  difference 
between  the  two  sets  of  data  shown  in  tables  3 and  4 is  that  time  lapse  values 
are  reduced  at  the  300  psi  level. 


(2)  Simple  Rectangle 


In  the  same  manner  as  the  previous  subsection,  general  descriptions 
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of  pressure  loading  around  the  rectangular  structure  are  presented  in  figures 
15  and  16.  These  figures  show  pressure  variations  as  a function  of  position 
for  specified  intervals  in  time.  Somewhat  different  behavior  characteristics 
are  evident  in  these  figures  than  were  found  in  figures  11  and  12  for  the  sim- 
ple arch  structure.  However,  the  five  basic  components  of  the  loading  phe- 
nomena-overpressure , reflections,  uniform  drag  effects,  vortex  behavior,  and 
rarefaction  effects  are  still  present. 

The  arrival  of  the  incident  shock  wave  results  in  high  reflected 
pressure''  on  the  front  face  of  the  structure  which  is  oriented  normal  to  the 
directio,  of  shock  propagation.  Unlike  the  arch  structure,  these  reflected 
pressures  are  maintained  for  a finite  period  of  time  related  to  the  location 
of  the  measurement  and  its  distance  from  the  top  edge  of  the  structure  face 
which  is  the  origin  of  rarefaction  waves.  Once  these  rarefaction  waves  arrive 
at  a given  point  the  reflected  pressures  begin  to  dissipate,  eventually  arriving 
at  a quasi-steady  pressure  level  composed  of  the  overpressure  and  a positive 
drag  component. 

On  the  top  surface  of  the  structure  each  location  is  initially 
loaded  by  the  incident  overpressure  level,  but  at  the  leading  edge  of  the  sur- 
face a vortex  is  generated  which  immediately  begins  to  erode  the  overpressure. 
After  a period  of  time  the  overpressure  is  completely  eliminated,  actually  re- 
sulting in  a slightly  negative  pressure  level  near  the  leading  edge.  At  posi- 
tions located  farther  downstream  from  the  leading  edge,  the  incident  pressure 
level  is  maintained  relatively  constant  for  a period  of  time  which  increases 
with  distance  from  the  leading  edge.  After  this  time  the  overpressure  level  is 
reduced  by  a negative  drag  pressure.  This  observation  of  increasing  time 
period  of  constant  overpressure  with  distance  from  the  leading  edge  suggests 
the  propagation  of  some  disturbance  from  the  leading  edge  with  a propagation 
velocity  less  than  the  incident  wave  speed.  Since  the  constant  overpressure 
level  is  generally  equal  to  the  incident  shock  wave  overpressure,  an  additional 
suggestion  is  that  the  incident  wave  initially  generates  a uniform  flow  with  a 
velocity  such  that  the  structure  surface  realizes  negligible  drag  effect.  At 
some  later  time  the  upstream  disturbance  arrives  and  alters  the  flow  resulting 
in  a negative  drag  pressure.  This  disturbance  may  well  be  the  expansion  of  the 
vortex  generated  at  the  leading  edge  of  the  top  surface.  The  accuracy  of  this 
statement  can  be  determined  only  by  experimental  optical  techniques.  The  data 
then  indicate  that  once  the  disturbance  arrives  at  a given  point  the  pressure 


POSITION 


POSITION 


100  Psi  Overpressure  Load  Distribution  - Rectangular  Structure 


Figure  15 


m«  erf  ur 

051  Q9T  I^T  111  21T 

I0SITIM 

031 

ii*  lit  e>r 

0)1  Ml  Ml  Ilf  Ilf 
WI1IH 

r 

600  r 

■ H ^uSlC 

- — — — 137  5 jiiec 

POSITION 





POSITION 


275  ^mc 
287.5^mc 
500  ^$»c 


237  5 n$*c 
2)0 

26?  3|UK 


05T  09T  I3T  I7T  2IT 
POSITION 


POSITION 


300  Psi  Overpressure  Load  Distribution 
Rectangular  Structure  (Continued) 


AFWL-TR-75-290 

loading  decreases  to  a relatively  stable  level. 

The  two  measurements  on  the  rear  of  the  structure  show  that  the 
initial  loading  wave  is  significantly  reduced  beneath  the  incident  overpressure 
level.  This  reduction  again  results  from  the  filling  process  of  the  incident 
wave  and  from  a strong  vortex  generated  at  the  top  of  the  rear  surface  by  the 
deflection  of  the  shocked  gas  around  the  corner.  At  later  times  the  loading 
increases  abruptly  as  a result  of  reflections  from  the  structure  support  plates. 
A similar  effect  would  be  seen  for  a structure  lying  on  a flat  surface  or  a 
structure  in  free  stream  conditions.  For  the  free  stream  condition,  a reflec- 
tion would  be  generated  when  the  incident  waves  from  the  top  and  bottom  collide. 
Even  with  the  reflections,  the  rear  surface  measurements  show  a pressure  level 
less  than  the  free  stream  overpressure  as  a result  of  the  negative  drag  pres- 
sure. 

Provided  in  tables  7 and  8 are  the  measured  peak  initial  pressures, 
steady  state  pressures,  and  time  lapse  values  as  functions  of  position  and  inci- 
dent pressure  level.  Calculated  incident  pressure  levels  derived  from  measured 
shock  front  velocities  are  also  shown.  These  tables  include  steady  state  pres- 
sure measurements  for  two  rear  face  locations.  However,  caution  should  be  used 
since  a stable  pressure  level  was  not  clearly  identified  on  the  rear  face  loca- 
tions, and  further,  the  recorded  data  at  later  times  in  this  area  include  re- 
flected pressure  effects  from  the  structural  model  support  plates. 

The  peak  initial  pressure  values  listed  for  front  face  gage  loca- 
cations  are  actually  peak  reflected  pressure  levels.  All  other  locations  indi- 
cate initial  wave  passage  over  each  gage  location. 

Tables  9 and  10  list  the  initial  peak  pressure  ratios  and  drag 
coefficients  calculated  from  the  measured  data  in  tables  7 and  8.  Average 
values  of  these  pressure  ratios  and  drag  coefficients  are  presented  graphically 
in  figures  17  and  18  for  each  incident  overpressure  level.  The  increase  in 
pressure  ratios  on  the  front  face  gage  locations  (reflection  factors)  with 
incident  overpressure  level  is  seen  in  figure  17  while  the  top  surface  shows 
the  essentially  undisturbed  transit  of  the  incident  shock  front.  The  rear 
face  locations  indicate  little  dependence  on  incident  pressure  level. 

Figure  18  suggests  that  an  increase  in  incident  pressure  level 
within  the  range  investigated  results  in  a slightly  lower  drag  coefficient. 

The  top  and  rear  surfaces  demonstrate  a more  pronounced  effect  with  lower  drag 
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Table  7 

RECTANGULAR  STRUCTURE  MEASUREMENTS- 100  PSI  OVERPRESSURE 
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Table  9 

RECTANGULAR  STRUCTURE  PRESSURE  RATIOS  AND  DRAG 
COEFF I C I ENTS- 1 00  PSI  OVERPRESSURE 


POSITION 

INCIDENT 

PRESSURE 

(PSI) 

INITIAL  PEAK 
PRESSURE  RATIO 

11  F 

94 

103 

5.64 

5.29 

09  F 

W~ 

103 

5.29 

5.10 

W~ 

5.85 

07  F 

103 

5.58 

94 

5.28 

05  F 

103 

5.49 

03  F 

94 

5.66 

03  T 

99 

1.09 

05  T 

99 

1.06 

09  T 

99 

1.05 

11  T 

99 

1.01 

13  T 

97 

T704- 

99 

0.98 

15  T 

97 

1.10 

99 

0.94 

17  T 

97 

1.01 

99 

1.11 

19  T 

99 

0.92 

21  T 

97 

99 

1.06 

0.99 

03  B 

no 

0.16 

05  B 

no 

0.16 

1.78 

1.70 

-0.81 

-0.73 

-0.68 

-0.54 

^039 

-0.63 

-0.49 

-0.54 

“3745 

-0.33 

-0.38 

“37R 

-0.48 

-0.47 

-0.45 
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Figure  17.  Initial  Peak  Pressure  Ratios  - Rectangular  Structure 
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coefficients  determined  for  the  100  psi  incident  overpressure  level. 

The  time  lapse  values  in  tables  7 and  8 show  increasing  values  on 
the  front  face  with  distance  from  the  top  edge.  This  is  to  be  expected  since 
the  duration  of  the  rarefaction  wave  increases  with  propagation  distance.  The 
increased  duration  results  from  the  faster  velocity  of  the  rarefaction  wave 
front  compared  to  the  propagation  velocity  of  the  tail.  On  the  top  surface  lo- 
cations the  time  lapse  values  generally  increase  with  distance  from  the  front 
edge  and  then  reach  a relatively  constant  value  on  the  latter  half  of  the  struc- 
ture surface.  This  trend  must  be  related  to  the  upstream  disturbance,  but  the 
nature  of  the  relationship  was  not  determined.  The  time  lapse  values  for  the 
rear  fare  measurvfnt nts  indicate  Irweesing  va'ftfot  with  JKtaroc  frtm  tt»  SUppBHi 
plates  as  one  might  expect  because  the  steady  state  pressures  are  pressure 
buildups  resulting  from  support  plate  reflection  and  drag  effects. 

The  pressure-time  histories  recorded  on  the  rectangular  structure 
front  face  were  of  a form  similar  to  that  shown  in  figure  19  where  three  specif- 
ic time  parameters  are  indicated. 


Figure  19.  Typical  Front  Face  Loading  - Rectangular  Structure 

The  value  Ti  represents  the  rise  time  to  peak  reflected  pressure 
Ppg  which  is  equivalent  to  the  rise  time  of  the  incident  overpressure.  The  sec- 
ond time  parameter,  T2,  represents  the  time  of  arrival  of  the  head  of  the  rare- 
faction wave  at  the  location  of  interest.  The  third  parameter,  T3,  represents 
the  time  at  which  the  tail  of  the  rarefaction  wave  passes  through  the  point  of 
interest  resulting  in  a steady  state  pressure  equal  to  the  free  stream  pressure 
and  a drag  component  of  loading.  The  value  of  T3  - T2  will  be  recognized  as  the 
time  lapse  value  discussed  earlier.  The  peak  reflected  pressure  duration,  T2  - 
Ti,  will  depend  upon  the  distance  of  the  point  of  interest  from  the  origin  of 
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the  rarefaction  wave  and  its  wave  front  velocity.  The  relief  time,  T3  - 12, 
will  then  depend  only  on  the  duration  of  the  rarefaction  wave.  Since  the  dura- 
tion of  the  rarefaction  wave  will  increase  with  propagation  distance,  relief 
time  should  be  found  to  increase  with  distance  from  the  origin  of  the  rarefac- 
tion wave.  For  the  rectangular  structure,  this  origin  is  in  the  vicinity  of  the 
top  edge.  With  these  considerations  in  mind  the  peak  reflected  pressure  dura- 
tion (Tj  - TJ  and  relief  ti.nts  [J,  - at  each  location  on  the  front  fate 
were  determined  and  plotted  as  shown  in  figure  20  for  each  of  the  incident  over- 
pressure levels.  These  values  were  rounded  off  to  the  nearest  five  microseconds. 
As  can  be  seen,  a somewhat  linear  relationship  exists  between  the  time  differ- 
ences and  the  distance  from  the  top  of  the  structure  front  face.  Additionally, 
the  curves  show  increasing  values  with  distance  consistent  with  the  previous 
discussion. 


(a)  100  psi  Overpressure 


O 

<L> 

l/l 


CVJ 


CO 


o 

Q) 

U) 


00 


(b)  300  psi  Overpressure 

Figure  20.  Peak  Reflected  Pressure  Durations  and  Relief  Times 
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Extension  of  the  curves  in  figure  20  into  the  left  vertical  axis 
which  is  the  location  of  the  top  edge  of  the  structure  face,  results  in  initial 
values  being  indicated  for  this  position.  The  initial  value  of  T2  - T15  if  it 
exists,  is  probably  related  to  the  expansion  of  reflected  pressure  into  the 
adjacent  region  near  the  edge  during  reflected  pressure  buildup  (T  < Tj).  The 
primary  rarefaction  wave  then  originates  at  some  point  away  from  the  edge  de- 
fined by  the  zone  of  expansion  at  T > Tj  when  increasing  reflected  pressures  no 
longer  occur  and  no  longer  prevent  the  propagation  of  the  rarefaction  wave  into 
the  zone  of  expansion.  The  initial  value  of  T2  - T i then  describes  the  transit 
time  of  the  rarefaction  wave  from  its  true  origin  to  the  edge  of  the  structure 
face. 

If  the  preceding  hypothesized  behavior  with  regard  to  T2  - Ti  is 
correct,  then  the  initial  value  will  be  dependent  only  on  the  shock  strength  of 
the  incident  overpressure  while  the  slope  of  the  T2  - Jx  curve  will  be  a func- 
tion of  the  rarefaction  wavefront  velocity  and  the  distance  of  the  point  of 
interest  from  the  edge  of  the  front  face.  An  estimate  of  the  rarefaction  wave- 
front  velocity  (CR)  may  be  obtained  from  the  following  expression  (ref.  3). 


CR  ' Co 


/Pso+Po\  /Pso+Po 

I po  M p 


PR0+Po \/PR0+P 


P +P 
so  0 


p +p 
SO  0 


<L+f^  1 

p T 1 


1 

' so  o' 


Using  average  values  of  peak  reflected  pressure  (PRQ)  across  the  front  face  and 
the  measured  incident  overpressure , rarefaction  wave  velocities  of  2320  ft/sec 
and  3720  ft/sec  were  determined  for  each  of  the  incident  overpressure  levels. 
These  values  compare  favorably  to  2430  ft/sec  and  3970  ft/sec  determined  from 
the  T2  - Ti  slopes  in  figures  20a  and  b respectively.  The  relationship  between 
T2  - T i and  position  may  therefore  be  expressed  as  follows: 


T2  - Tj  = A + 


where  A represents  the  initial  time  valffe  at  the  edge  location  and  S represents 
the  distance  of  the  point  of  interest  from  the  edge  of  the  structure  face.  For 
realistic  structure  sizes  the  value  of  A becomes  inconsequential  so  that  the 
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peak  reflected  pressure  duration  could  be  tan  , equal  to  f— . 

lR 

The  T3  - T2  versus  position  curves  in  figure  20  for  each  of  the 
overpressures  show  the  increase  in  rarefaction  wave  duration  with  propagation 
distance.  These  durations  may  also  be  expressed  as  an  initial  value  plus  an 
expression  in  terms  of  the  propagation  distance  and  rarefaction  wave  speed. 

The  expressions  result  in  values  of  1.41  S/CR  and  1.34  S/CR  for  each  of  the 
overpressure  levels  investigated.  Since  the  initial  values  of  the  T3  - T2 
curves  at  the  edge  of  the  structure  face  are  relatively  large,  these  values 
cannot  be  ignored  for  structural  dimensions  considered  in  this  investigation. 
However,  for  much  larger  dimensions,  the  importance  again  becomes  minimal  with 
respect  to  the  slope  equation,  assuming  of  course  that  the  slope  equation  re- 
mains valid  for  larger  dimensions. 

b.  Extended  Roof  Models 

The  objective  of  the  extended  roof  tests  was  to  determine  whether  the 
front  and  top  surface  loadings  would  be  independent  of  the  downstream  exten- 
sions. The  approach  taken  to  satisfy  this  objective  was  to  measure  surface 
loadings  at  the  70,  80,  and  90  degree  locations  on  the  arch  structure  and  at 
the  13T,  15T,  17T,  19T,  and  21T  locations  on  the  rectangular  structure;  both 
structures  had  the  roof  extensions  in  place.  These  measurements  were  then  com- 
pared to  earlier  measurements  taken  without  the  extensions,  but  at  similar  gage 
locations.  These  comparisons  are  presented  in  the  following  subsections. 

(1)  Extended  Arch 

Comparisons  of  the  pressure-time  histories  recorded  for  the  simple 
arch  structure  (SA)  and  the  extended  arch  structure  (EA)  are  presented  in 
figures  21  and  22.  Figures  21  a,  b,  and  c show  the  100  psi  incident  overpres- 
sure measurements  at  the  70,  80,  and  90  degree  locations  respectively.  The  com- 
parison presented  for  the  70  degree  location  supports  the  contention  of  no 
effect  due  to  downstream  extensions.  However,  the  80  and  90  degree  locations 
indicate  otherwise,  showing  a more  stable  waveform  for  the  extended  models. 
Obviously,  the  initial  difference  in  magnitudes  shown  in  figures  21  b and  21  c 
can  result  only  from  differences  in  gage  sensitivity,  but  this  sensitivity 
difference  would  result  in  consistent  magnitude  differences  throughout  the 
period  of  recording  and  the  waveforms  would  be  similar.  The  contradiction 
shown  by  the  dissimilarity  in  waveforms  suggests  that  the  downstream  extensions 
do  affect  the  pressure  loadings  in  areas  of  proximity  to  the  extension. 
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This  effect  would  seem  to  be  a geometry  change  relationship  which  will  tend  to 
alter  the  flow  around  the  structure  and  the  resulting  dynamic  pressures. 

Arch  structure  loading  comparisons  for  300  psi  incident  overpres- 
sures are  presented  in  figures  22  a,  b,  and  c.  The  70  degree  locations  exhibit 
no  pronounced  differences  in  either  waveform  or  magnitude  whereas  the  80  and 
90  degree  locations  again  seem  to  realize  the  effects  of  the  downstream  exten- 
sion. Also  quite  evident  in  these  figures  are  later  effects  of  combustion 
product  arrivals  for  the  extended  models. 

(2)  Extended  Rectangle 

Simple  rectangle  (SR)  and  extended  rectangle  (ER)  pressure-time 
histories  are  compared  in  figures  23  and  24  for  the  100  psi  and  300  psi  over- 
pressure levels.  The  obvious  difference  shown  in  the  comparisons  is  for  the 
extended  model  data  to  display  somewhat  higher  pressure  values  after  arrival 
of  the  initial  loading  wave.  The  higher  values  are  in  most  cases  then  sus- 
tained throughout  the  recording  time  (prior  to  combustion  products  arrival). 
This  result  suggests  that  the  extension  produces  a slight  alteration  in  the 
flow  around  the  rectangular  structure  which  results  in  a higher  total  pressure 
acting  on  the  top  surface  at  the  locations  considered.  The  comparison  shown  in 
figure  23  e shows  that  the  location  near  the  top  edge  (100-SR-21T)  realizes  a 
significant  effect  from  the  deflection  of  the  shock  into  the  rear  face  region. 
However,  the  similar  data  for  the  300  psi  level  shown  in  figure  24  e show  a 
negligible  effect.  These  two  observations  suggest  that  the  300  psi  flow  veloc- 
ity exceeds  the  propagation  velocity  of  any  major  disturbance  from  the  down- 
stream geometry  change.  The  deviation  in  waveform  shown  in  figure  24  d for  the 
extended  model  19T  location  does  not  seem  to  be  a result  of  the  model  extension 
since  a similar  disturbance  was  not  noted  at  the  21T  location.  This  same 
reasoning  indicates  that  the  deviation  does  not  result  from  an  upstream  distur- 
bance cither,  based  on  the  traces  shown  in  figures  24  a,  b,  and  c.  Evidently 
the  waveform  change  results  either  from  a localized  disturbance  or  from  tran- 
sient malfunctioning  of  the  gage  or  recording  system. 

c.  Wedge  Extensions  (Simulated  Ground  Planes) 

The  simulated  ground  plane  tests  were  designed  to  collect  data  which 
would  describe  the  behavior  of  reflected  waves  propagating  upstream  'from  the 
test  structures.  These  data  consisted  of  pressure-time  measurements  at  the 
stagnation  points  (midheight)  on  the  test  structures  and  at  one  half  and  one 
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structure  height  in  front  of  the  arch  and  rectangle.  It  was  originally  antici- 
pated that  the  pressure-time  recordings  would  additionally  yield  information 
concerning  rarefaction  wave  behavior,  steady  state  pressure,  and  relief  times 
to  steady  state  loading.  The  requirements  for  this  information  were  based  on 
a preconception  that  the  pressure-time  history  at  a point  away  from  the  struc- 
ture would  be  of  the  form  shown  in  figure  25.  This  waveform  would  describe 
the  incident  wave  passage  followed  by  the  reflected  wave  from  the  structure. 

It  was  then  assumed  the  reflected  wave  would  decrease  uniformly  to  a steady 
state  level  composed  of  free  stream  overpressure  and  an  influence  from  the 
dynamic  pressure  effects  on  the  front  face  of  the  structure.  As  the  experimen- 
tal data  in  figures  26  and  27  demonstrate,  this  preconception  was  somewhat  in 
error.  These  figures  show  pressure-time  measurements  at  the  stagnation  point 
and  at  one  half  and  one  structure  height  in  front  of  the  structures  for  each 
of  the  incident  overpressure  levels.  Instead  of  a simple  reduction  to  steady 
state  conditions  after  arrival  of  the  reflected  wave,  the  data  show  arrivals 
of  secondary  loading  effects  from  the  structure. 
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The  100  psi  arch  data  in  figure  26  show  a marked  attenuation  of  the 
reflected  wave  at  the  0.5R  position  followed  by  a gradual  rise  in  pressure 
which  may  be  attributed  to  the  dynamic  pressure  effects  acting  on  the  structure 
face  and  secondary  reflections  froir.  the  region  near  the  intersection  of  the 
simulated  ground  plane  and  the  structure  face.  Such  reflections  have  been 
identified  in  Eulerian  code  calculations  of  shocks  propagating  into  closed 
end  tunnels  where  the  initial  reflected  wave  from  the  end  face  was  found  to 
reflect  off  the  adjacent  tunnel  walls.  The  1.0R  location  shows  negligible 
attenuation  of  the  reflected  wave  from  the  0.5R  location.  The  initial  re- 
flected wave  magnitude  at  the  1.0R  location  is  then  sustained  relatively  con- 
stant for  the  remaining  recording  time.  This  sustained  level  is  probably  the 
result  of  the  dynamic  pressure  effects  generated  by  the  structure. 

The  data  in  figure  26  b for  the  300  psi  arch  show  similar  behavior. 
Again,  a significant  attenuation  of  the  reflected  wave  occurs  between  the 
structure  face  and  the  0.5R  location.  Between  0.5R  and  1. OR,  only  a small  re- 
duction in  reflected  pressure  is  seen.  Pressure  levels  at  both  locations  then 
remain  essentially  constant.  At  the  300  psi  level,  the  0.5R  location  does  not 
display  the  gradual  increase  in  pressure  noted  for  the  100  psi  level.  This 
observation  suggests  that  the  300  psi  incident  pressure  level  is  less  suscep- 
tible to  high  stccndary  reflection!,  than  the  100  psi  level.  The  fact  that  tuth 
positions  register  fairly  constant  pressure  levels  may  indicate  that  the  sec- 
ondary reflections  are  within  the  pressure  levels  measured  for  the  initial 
attenuated  reflection  and  the  later  dynamic  pressure  effects  from  the  structure. 

The  data  collected  for  the  100  psi  rectangular  structure  (figure  27  a) 
show  lower  attenuation  rates  and  more  pronounced  secondary  loads  with  late  time 
pressures  converging  on  the  steady  state  pressure  level  measured  at  the  stagna- 
tion point.  The  lower  attenuation  might  be  expected  since  the  high  initial 
reflected  pressure  is  generated  over  a larger  region  and,  further , these  peak 
reflected  values  are  sustained  for  a lonjpr  period  of  tiW3  On  the  fac*>  uf  the 
structure.  The  convergence  of  both  curves  to  the  steady  state  level  measured 
at  the  stagnation  point  suggests  that  a constant  dynamic  pressure  effect  acts 
over  a finite  region  in  front  of  the  rectangular  structure.  For  the  100  psi 
level  this  region  exceeds  the  structure  height. 

For  the  arch  structure  at  the  100  psi  level  this  region  is  approximate- 
ly one  half  the  structure  height  while  at  300  psi  this  region  is  noticeably 
less  than  one  half  the  structure  height. 


AFWL-TR- 75-290 

The  300  psi , rectangular  structure  data  shown  in  figure  27  b exhibit 
a greater  attenuation  of  the  initial  reflected  wave  between  the  structure  and 
the  0.5H  location  than  was  shown  for  the  100  psi  case.  Between  0.5  and  1.0H 
the  change  in  reflected  pressure  level  is  small  followed  by  a rapid  decrease 
to  a relatively  stable  pressure  level.  After  a period  of  time,  the  pressure 
again  increases.  This  increase  is  too  advanced  in  time  to  be  a result  of 
combustion  products  yet  it  seems  too  late  to  be  a result  of  secondary  reflec- 
tions or  dynamic  pressure  effects.  The  true  source  of  the  increase  was  not 
identified.  The  region  of  constant  dynamic  pressure  effects  from  the  struc- 
ture extends  beyond  the  0.5H  location,  but  does  not  reach  the  1.0H  location. 
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SECTION  V 

EVALUATION  OF  CURRENT  PREDICTION  TECHNIQUES 


1 . INTRODUCTION 

The  primary  objective  of  this  investigation  was  to  determine  the  adequacy 
of  existing  methods  to  predict  loadings  on  aboveground  structures  subjected  to 
high  incident  overpressures.  In  this  report  the  investigation  was  limited  to 
those  prediction  techniques  of  c simplified  nature.  These  simplified  proce- 
dures allow  a given  overpressure -time  history  to  be  transformed  into  surface 
loads  on  a structure  through  the  use  of  such  parameters  as  reflection  factors 
and  drag  coefficients  and  also  with  the  use  of  shock  wave  characteristics  and 
structure  geometry  considerations.  The  Air  Force  Manual  for  the  Design  anH 
Analysis  of  Hardened  Structures  (ref.  3)  was  used  as  a source  reference  for 
the  prediction  techniques  evaluated.  This  document  represents  the  most  current 
collection  of  recommended  procedures  for  preliminary  design  or  analysis  of 
hardened  strategic  structures. 

The  format  followed  in  evaluating  these  techniques  was  first  to  compare  the 
basic  transformation  parameters  proposed  by  the  previous  document  with  those 
determined  in  this  effort.  The  actual  structural  loads  measured  during  the 
investigation  were  then  compared  with  those  resulting  from  the  prediction 
techniques . 

2.  REFLECTION  FACTORS 

a.  Stagnation  Probe  Tests 

Reproduced  in  figure  28  is  the  experimentally  determined  reflection 
factor  curve  from  the  stagnation  probe  test  data.  Superimposed  on  this  figure 
is  the  normal  reflection  factor  curve  recommended  by  reference  3.  The  compari- 
son shows  higher  values  for  the  experimental  data.  These  higher  values  should 
be  expected  since  the  experiments  were  conducted  under  acmospheric  pressures 
of  approximately  12  psi  whereas  the  recommended  curve  was  determined  for  sea 

level  air.  The  comparison  can  therefore  lend  only  additional  support  to  the 
recommended  curve. 
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Figure  28.  Variation  of  Reflection  Factor  with  Incident 
Overpressure  - Stagnation  Probe  Tests 


b.  Simple  Arch  Structure 

For  the  overpressures  of  interest  in  this  investigation  the  "Design 
Manual"  does  not  clearly  specify  recommended  reflection  factors  as  functions  of 
angle  of  incidence  on  an  arch  structure.  The  "Manual"  does,  however,  provide 
reflection  factors  for  shock  waves  impinging  on  inclined  planar  surfaces.  Since 

I a user  of  this  document  will  likely  attempt  to  use  these  values  in  any  analysis 

requiring  point  load  determinations,  the  comparisons  in  figure  29  are  presented 
for  evaluation.  T.he  term  "recommended"  is  used  with  a great  degree  of  license 
in  these  comparisons. 


The  recommended  and  experimental  curves  for  the  100  psi  level  as  a 
whole  demonstrate  good  agreement  between  reflection  factors.  The  experimental 
curve  differs  primarily  by  showing  a greater  variability  in  reflection  factors 
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Th  ea  r 1 , . 2er°  a"d  the  Mtre,ne  a"9le  for  Mach  reflection 

J\  local,0"  of  the  extreme  angle  is  not  clearly  evident  on  the  experi 

curve  as  a result  of  the  poor  resolution  obtained  with  the  relatively 

coarse  gage  spacing.  The  true  location  of  the  extreme  angle  probably  lies 

either  at  the  40  degree  location  or  at  a slightly  higher  value  since  the  30 

d 40  degree  values  represent  a positive  change  in  slope.  If  the  extreme 

angle  was  at  some  location  between  30  and  40  degrees,  a pronounced  negative 
slope  change  would  be  indicated. 


values  . ' , , f ,ndde"Ce  ”i9her  tha"  the  eXtrere  a"9le>  ^ experimenta 

Slues  show  Slightly  decreased  reflection  factors.  This  observation  has  even 

greater  impact  remembering  that  the  experimental  values  were  determined  under 

ambi  nt  pressures  of  approximately  12  psi  while  the  recommended  values  are  fo, 

« taTdit?"5  Whkh  W0U,d  y,eld  W factors  for  the 
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The  comparisons  for  the  300  psi  level  indicate  two  areas  of  disagree- 
ment between  the  data.  First,  the  experimental  values  are  lower  in  magnitude 
within  the  region  between  zero  angle  of  incidence  and  the  extreme  angle,  again 
remembering  the  effect  of  lower  ambient  pressure.  Second,  the  extreme  angle 
is  clearly  located  at  an  angle  in  excess  of  40  degrees  resulting  in  higher 
reflection  factors  encompassing  a greater  portion  of  the  structure. 

c.  Simple  Rectangular  Structure 

Because  of  the  nature  of  its  geometry  and  orientation  to  the  shock  wave, 
the  simple  rectangle  will  have  an  essentially  constant  reflection  factor  over 
its  entire  front  face  for  any  given  overpressure  level.  Using  composite  pres- 
sure-time histories  which  yield  the  average  pressure  on  the  front  face,  reflec- 
tion factor  values  of  4.8  and  6.4  were  determined  for  the  100  and  300  psi  over- 
pressure levels.  These  values  compare  favorably  with  reflection  factors  of 
4.9  and  6.6  recommended  by  the  "Design  Manual." 

3.  DRAG  COEFFICIENTS 

a.  Stagnation  Probe  Tests 

A comparison  of  the  recommended  drag  coefficient  curve  for  a disc- 
shaped object  and  the  experimentally  determined  curve  is  presented  in  figure 
30.  The  experimental  curve  shows  a more  definite  incre,.'.e  in  drag  coefficient 
with  increasing  overpressure  level.  This  larger  increase  may  not  be  a com- 
pletely accurate  representation  because  of  the  small  number  of  test  conducted 
at  the  higher  overpressure  levels.  These  values  are  all  within  the  range  of 
scatter  experienced  at  the  lower  pressure  levels  (see  figure  10).  Only  addi- 
tional testing  will  determine  truly  representative  values  at  these  higher  pres- 
sure levels.  The  larger  number  of  lower  pressure  measurements  does  support  the 
conclusion  that  the  recommended  curve  (for  lower  pressures)  overestimates  drag 
coefficient  values  for  these  test  conditions  and  the  manner  in  which  the 
dynamic  pressures  were  determined. 

b.  Simple  Arch  Structure 

The  "Design  Manual"  does  not  provide  recommended  drag  coefficients  as  a 
function  of  angle  of  incidence  for  the  overpressure  levels  investigated  in  this 
report.  The  "Manual"  does  however  suggest  that  for  simplified  loading  conditions 
the  total  drag  coefficient  for  an  arch  structure  at  high  overpressures  may  be 
taken  equal  to  one.  This  value  compares  fairly  well  with  values  of  1.25  and 
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1.15  determined  from  composite  records  which  were  constructed  from  the  hori 
tal  pressure  components  acting  at  each  of  the  gage  locations.  * 


EXPERIMENTAL 

RECOMMENDED  (Ref.  3) 


14C  190  240 

INCIDENT  OVERPRESSURE  (PSI) 

Figure  30.  Variation  of  Drag  Coefficient  with  Incident  Over- 
pressure - Stagnation  Probe  Tests 

c.  Simple  Rectangular  Structure 

Surface  drag  coefficients  are  presented  in  the  "Design  Manual"  for  a 
rectangular  structure,  but  only  for  average  loads  on  an  entire  surface  such 
as  the  front,  roof,  or  rear  face.  For  the  roof  and  rear  surfaces  these  recom- 
mendations are  further  restricted  to  incident  overpressures  less  than  130  psi. 
In  order  to  derive  similar  drag  coefficient  values,  it  was  necessary  to  use 
composite  records  constructed  from  measurements  common  to  each  of  the  different 
structure  surfaces.  These  composites  are  presented  in  subsection  4 where 
total  loads  on  the  structures  will  be  presented. 

Recommended  front  face  drag  coefficients  of  1.25  and  1.45  are  specified 
for  the  100  and  300  psi  overpressure  levels.  The  experimental  data  yielded 
values  of  1.3  and  1.5  for  each  of  these  pressure  levels  which  resulted  in  ex- 
cellent correlation  between  the  experimental  and  recommended  values. 

A value  of  -0.2  is  recommended  for  a drag  coefficient  on  the  roof  of 
the  rectangular  structure  at  an  overpressure  level  of  50  to  130  psi.  The  ex- 
perimental data,  on  the  other  hand,  indicate  a drag  coefficient  of  approxi- 
mately -0.6  at  the  100  psi  level.  At  the  300  psi  level  the  drag  coefficient 
changes  to  a value  of  -0.3. 
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(b)  300  psi  Overpressure 

Figure  32.  Predicted  Average  Surface  Loadings-Arch  Structure 
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Experimental  drag  coefficients  of  -0.5  and  -0.1  were  determined  for 
the  100  and*  300  psi  overpressure  levels  on  the  rear  face.  However,  these 
values  are  not  truly  representative  of  the  entire  rear  face  since  only  two 
pressure  records  were  obtained  at  locations  on  the  upper  portion  of  the  rear 
face  (03B  and  05B).  The  "Design  Manual"  recommends  a value  of  -0.2  for  over- 
pressure levels  between  50  and  130  psi. 


It  would  seem  that  roof  and  rear  face  drag  coefficients  should  be 
dependent  upon  surface  size  when  specified  as  averages  for  an  entire  face. 

This  statement  is  made  in  view  of  the  occurrence  of  vortex  behavior  and  pos- 
sible flow  separation  at  the  leading  edges  of  both  the  roof  and  rear  faces. 

The  accompanying  reductions  in  pressure  will  therefore  have  a greater  effect  as 
the  surface  dimension  in  the  direction  of  flow  around  the  structure  decreases. 


4.  TOTAL  STRUCTURAL  LOADS 


A common  technique  used  in  preliminary  design  or  analysis  is  to  assume 
that  surface  loads  on  a structure  are  uniform  along  any  given  face  of  the 
structure  such  that  time  is  the  only  variable  considered.  This  technique 
allows  estimates  to  be  made  of  overturning  moments,  translational  forces,  and 
foundation  reactions.  In  actuality,  these  estimates  are  often  used  for  final 
design  or  analysis  for  lack  of  more  precise  capabilities.  It  is  the  purpose 
of  this  section  to  compare  simplified  analysis  results  with  those  determined 
experimentally. 


a.  Simple  Arch  Structure 


Under  simple  loading  conditions  the  "Design  Manual"  suggests  an  analy- 
sis which  assumes  an  arch  structure  will  respond  in  two  primary  modes;  a 
breathing  mode  corresponding  to  a uniform  overpressure  and  a flexural  mode 
corresponding  to  a load  applied  radially  inward  on  the  upstream  face  and 
radially  outward  on  the  downstream  face  (see  fig.  31).  The  assumed  flexural 
loading  is  composed  of  two  components  which  consider  the  unsymmetrical  loads 
imparted  to  the  structure  as  the  shock  wave  passes  over  the  arch  and  a drag 
component  resulting  from  the  continuing  drag  loading  after  the  arch  has  been 
engulfed  by  the  shock  wave.  The  uniform  radial  pressure  is  assumed  to  increase 


linearly  from  zero  to  the  maximum  overpressure  level  in  a time  t given  by 


*r  = O - 0/TT ) r 
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Figure  31.  Assumed  Response  Modes  for  Arch  Structure 


where 


t = transit  time  of  the  shock  wave  across  the  structure 
<f>  = half  the  central  angle  of  the  arch 

Thereafter,  the  uniform  radial  pressure  is  assumed  to  vary  as  the  free  stream 
overpressure.  The  initial  component  of  the  flexural  loading  is  assumed  to 

(1/2  + cp/n)  P 


C ,q 
d o 


so  TS7fl+6_ 


increase  linearly  from  zero  to  a maximum  value  of 

at  a time  r/2  and  to  decay  linearly  to  zero  at  a time  (1  + 3<j>/Tr)r.  The  drag 
component  of  the  flexural  loading  increases  linearly  from  zero  to  a maximum 
value  of  (<f>/iT)cdqo  at  a time  (1  + 30/tt)t  and  then  remains  equal  to  (cp/iT )cdq (t) 
For  this  effort  the  peak  dynamic  pressure,  qQ,  and  time  variation  of  dynamic 
pressure,  q(t),  are  equivalent.  The  "Design  Manual"  suggests  that  for  high 
incident  overpressures,  a value  of  1 may  be  used  for  the  drag  coefficient, 
c 


“d’ 


Using  the  above  techniques,  time  variations  of  pressure  as  shown  in 
figure  32  were  determined  for  the  100  and  300  psi  incident  overpressure  levels 
The  pressure  levels  at  any  given  time  shown  in  these  figures  may  be  compared 
to  the  pressure  variations  at  similar  times  shown  earlier  in  figures  11  and  12 
where  pressures  as  a function  of  position  were  presented.  In  making  this  type 
of  comparison  it  is  evident  that  not  only  do  the  predicted  curves  oversimplify 
the  loading  conditions,  they  also  underestimate  the  average  loads  acting  on 
either  face  of  the  arch  at  any  given  time.  The  effect  of  these  differences  on 
foundation  reactions  may  be  determined  by  considering  the  restraints  necessary 
to  prevent  rigid  body  motion  excluding  inertial  effects.  The  arch  will  be 
assumed  rigid  and  to  have  pinned  connections  at  the  0 and  180  degree  points. 
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Presented  in  figure  33  are  the  average  vertical  and  horizontal  pres- 
sures acting  on  the  arch  and  the  computed  overturning  moment  about  the  down- 
stream connection  for  the  100  psi  overpressure  level.  During  the  true  diffrac- 
tion phase  the  predicted  horizontal  loading  underestimates  the  actual  by 
approximately  30  to  40  percent  which  is  equivalent  to  stating  that  the  required 
lateral  restraint  would  be  similarly  underestimated.  After  the  diffraction 
Phase  the  differences  between  actual  and  predicted  levels  are  small.  The  com- 
parison in  figure  33  b between  the  experimental  and  predicted  average  vertical 
pressures  shows  excellent  agreement.  However,  this  measurement  is  meaningless 
unless  it  is  used  with  the  calculated  moments  in  figure  33  c to  determine  the 
vertical  reactions  at  each  of  the  connection  points.  The  close  correlations 
between  the  experimentally  determined  moments  and  predicted  moments  and  the 
agreement  shown  for  the  average  vertical  loadings  suggest  that  the  prediction 

technique  is  satisfactory  for  determining  vertical  reactions  at  the  structure 
supports. 

Similar  pressure  and  moment  determinations  are  shown  in  figure  34  for 
the  300  psi  overpressure  level.  The  predicted  horizontal  pressures  again 
underestimate  the  measured  quantities  showing  greater  differences  during  early 
times  than  were  found  for  the  100  psi  level.  The  late  time  quasi-steady  pres- 
sure level  is  underestimated  by  approximately  15  percent.  Unlike  the  previous 
lower  pressure  level,  the  300  psi  data  show  dissimilar  comparisons  for  the 
average  vertical  pressures  (fig.  34  b).  These  differences  plug  the  major 
discrepancies  between  predicted  and  experimental  overturning  mounts  shown  in 
figure  34  c will  result  in  erroneous  levels  of  predicted  structure  support 
reactions.  After  the  maximum  predicted  levels  of  vertical  pressure  and  moment 
are  reached,  the  predicted  levels  result  in  an  average  underestimate  of  30 
percent  at  the  upstream  connection  and  an  average  overestimate  of  18  percent 
at  the  downstream  support  (vertical  reaction  loads  only).  These  differences 
P us  the  disagreement  fiund  between  predicted  and  measured  average  horizontal 
pressures  suggest  that  the  simplified  analysis  techniques  are  somewhat  unfavor- 
able  for  the  higher  300  psi  overpressure  level. 

Both  the  100  psi  and  300  psi  predicted  time  histories  of  horizontal 
pressure  resulted  in  large  underestimates  of  average  loading  during  the  dif- 
fraction phase.  A review  of  the  simplified  analysis  techniques  presented 
earlier  shows  that  these  underestimates  would  quite  naturally  result  since 
these  techniques  do  not  accurately  address  the  reflection  behavior  on  the  front 
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(c)  Overturning  Moment 

Figure  34.  Experimental  and  Predicted  Arch  Loadings  - 
300  Psi  Overpressure  (Continued) 

face.  The  maximum  average  pressure  on  the  front  face  for  any  given  overpres- 
sure is  limited  by  the  value  of  the  [1/2  + term  in  the  expression  for  the 
initial  flexural  loading.  For  the  cases  investigated  this  term  has  a value 
of  1.  Since  the  drag  component  of  flexural  loading  represents  a negligible 
quantity  during  the  time  of  peak  loading,  the  peak  average  reflected  pressure 
is  effectively  limited  to  2Psq  on  the  front  face.  Referring  back  to  figures 
11  and  12  the  data  indicate  that  the  average  peak  values  are  on  the  order  of 
3Psq  (300  psi  and  1000  psi)  for  each  of  the  overpressures  investigated.  The 
error  is  compounded  by  the  assumption  of  uniform  pressure  across  the  surface. 
The  higher  pressure  levels  occur  at  the  lower  angles  of  incidence  (less  than 
= 45°)  where  the  loading  is  predominant ly  horizontal.  The  uniform  assumption 
would  result  in  a lower  horizontal  conponent  even  if  the  average  pressure  level 
was  accurate.  In  a similar  sense  the  uniformity  assumption  also  will  affect 
the  values  determined  for  overturning  moment. 

The  later  time  quasi-steady  loading  showed  somewhat  reasonable  predic- 
tion values  in  terms  of  horizontal  loading.  However,  as  shown  in  figures 
33  c and  34  c,  the  predicted  overturning  moments  differ  appreciably  from  the 
measured  data.  Again,  this  difference  is  directly  related  to  the  assumption  of 
uniform  surface  loading  on  the  front  face.  In  all  cases  the  assumption  of 
uniform  pressure  on  the  rear  face  is  not  too  unreasonable  since  the  pressure 
levels  tend  to  vary  about  a median  level. 
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For  the  average  vertical  loads  across  the  entire  structure,  the 
prediction  techniques  result  in  uniform  pressure  equal  to  the  incident  over- 
pressure. This  result  is  caused  by  the  cancellation  of  vertical  components 
from  the  flexural  loadings  which  are  radial  and  in  on  the  front  face  and  radial 
and  out  on  the  rear  face.  In  real i ty,  this  effect  could  be  realized  during 
the  post  diffraction  loadings  and  under  simple  loading  conditions  only  if  the 
average  drag  coefficients  on  the  front  face  were  equal  in  magnitude  and  oppo- 
site in  sign  to  the  rear  face  value.  Figure  14  shows  this  not  to  be  the  case 
and  figures  33  b and  34  b show  the  result  at  late  times.  For  low  overpressure 
levels  the  assumed  behavior  may  be  more  accurate. 

The  predicted  pressure -time  histories  additionally  underestimate  the 
time  duration  of  the  diffraction  process  which  is  limited  by  the  [1  + 3<J)/tt] 
term.  A more  accurate  value  for  the  conditions  examined  is  [1  + 7<J)/Tr]. 


b.  Simple  Rectangular  Structure 

The  "Design  Manual"  techniques  for  estimating  average  loads  on  an 
aboveground  rectangular  structure  are  based  on  the  simplified  loading  behavior 
shown  in  figure  35.  The  average  load  on  the  front  face  is  assumed  to  rise  to 
a peak  reflected  pressure  (PRq)  in  a time  (TR)  equal  to  the  rise  time  of  the 
incident  overpressure.  The  value  of  PRq  is  determined  from  the  incident  over- 
pressure and  its  corresponding  reflection  factor  for  normal  reflection.  The 
peak  reflected  pressure  is  then  assumed  to  decay  linearly  to  a quasi -steady 
pressure  level  equal  to  the  overpressure  plus  the  drag  component  of  loading 
which  is  determined  from  the  free  stream  dynamic  pressure  and  front  face  drag 
coefficient.  This  quasi-steady  loading  is  assumed  to  decay  in  the  same  manner 
as  the  free  stream  overpressure  plus  resulting  drag  pressure.  The  decay  time 
of  the  reflected  pressure  is  assumed  equal  to  the  time  required  for  a rare- 
faction wavefront  to  travel  from  the  nearest  edge  of  the  front  face  to  the 
stagnation  point  and  then  return  to  the  edge  of  the  structure  (stagnation 
time).  The  average  roof  loading  on  the  structure  is  assumed  to  rise  linearly 
to  a peak  value  in  a time  equal  to  the  rise  time  of  the  blast  wave  plus  the 
transit  time  of  the  shock  wave  across  the  roof.  The  peak  value  is  defined  as 
being  equal  to  the  incident  overpressure  reduced  by  a drag  component  resulting 
from  i negative  drag  coefficient  and  the  peak  dynamic  pressure.  After  this 
time  the  average  pressure  is  taken  equal  to  the  current  overpressure  reduced 
by  the  negative  drag  pressure  determined  from  the  current  dynamic  pressure. 

The  "Design  Manual"  additionally  provides  a similar  technique  for  determining 
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the  average  loading  on  the  rear  face  of  the  rectangular  structure.  However, 
an  insufficient  number  of  measurements  were  conducted  on  this  face  to  justify 
a comparison. 


(a)  Front  Face 


(b)  Roof 


Figure  35.  Simplified  Loading  Behavior  - Rectangular  Structure 

Using  the  above  techniques,  average  pressures  were  determined  for  the 
front  and  roof  surfaces  of  the  rectangular  structure  and  are  compared  to  the 
experimentally  determined  averages  in  figures  36  and  37.  The  predicted  time 
histories  for  the  100  psi  level  shown  in  figure  36  were  determined  using  a 
reflection  factor  of  4.9,  a front  face  drag  coefficient  of  1.25',  and  a roof 
drag  coefcicient  of  -0.2.  As  the  comparisons  demonstrate,  the  prediction 
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Figure  36. 


Pr®dicted  Rectangular  Structure 
Loadings  - 100  Psi  Overpressure 
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techniques  yield  favorable  results  for  the  front  face  loading,  with  the  excep- 
tion of  stagnation  time,  and  unsatisfactory  results  for  the  roof  loading.  The 
experimental  data  show  that  for. the  front  face  more  appropriate  load  parameters 
would  be  a reflection  factor  of  4.8,  a drag  coefficient  of  1.3,  and  a stagna- 


tion time  of  4S/Cr.  The  increased  stagnation  time  is  probably  the  result  of 


the  simplified  procedure  not  taking  into  account  the  duration  of  the  rarefac- 
tion wave.  Excellent  agreement  is  shown  for  the  roof  loading  during  the  tran- 
sit time  of  the  shock  across  the  structure  top  surface.  However,  the  simpli- 
fied technique  significantly  overestimates  the  average  loading  occurring 
during  the  quasi-steady  loading  period.  As  discussed  earlier,  the  initial 
flow  across  the  top  surface  is  probably  parallel  to  that  surface  such  that  a 
drag  coefficient  of  -0.2  is  appropriate  for  determining  the  initial  peak  pres- 
sure level.  However,  once  the  flow  field  has  stabilized  following  the  diffrac- 
tion period  the  flow  is  no  longer  parallel  across  the  top  surface  of  the  struc- 
ture. This  resulting  nonparallel  flow  results  in  a quasi-steady  drag  coeffi- 
cient of  -0.6.  Assuming  that  the  relationship  is  scalable,  the  experimental 
findings  suggest  that  for  structures  with  similar  height  (stagnation  distance) 
to  length  ratio  (S/L  = 2)  the  average  roof  loading  may  be  described  by  a 
linear  rise  to  a peak  pressure  defined  by  the  overpressure  reduced  by  a nega- 
tive drag  component  resulting  from  a drag  coefficient  of  -0.2  in  a time  equal 
to  the  transit  time  across  the  roof  surface.  The  quasi -steady  loading  may  then 
be  taken  equal  to  the  current  overpressure  reduced  by  a negative  drag  component 
determined  by  a drag  coefficient  of  -0.6.  The  decrease  from  peak  pressure 
level  to  quasi-steady  level  may  be  assumed  linear  and  to  occur  over  a period 
of  time  equal  to  twice  the  shock  wave  transit  time  across  the  surface. 


The  predicted  front  and  top  surface  average  pressure  loadings  are  com- 
pared in  figure  37  to  the  measured  time  histories  for  the  300  psi  overpressure 
level.  A reflection  factor  of  6.5,  front  face  drag  coefficient  of  1.45  and  a 
rear  face  drag  coefficient  of  -0,2  were  used  in  the  predictions.  The  front 
face  loading  comparison  demonstrates  excellent  agreement  between  predicted  and 
measured  values.  The  only  discrepancy  shown  is  that  the  predicted  curve 
indicates  a shorter  clearing  time  than  was  actually  measured.  The  actual 
stagnation  time  was  approximately  3.5  S/CR.  However,  to  maintain  the  total 
impulse  delivered  to  the  front  face  during  the  diffraction  period,  a linear 
decay  to  quasi-steady  pressure  level  in  a time  equal  to  3 S/CR  is  more  appro- 
priate for  the  300  psi  overpressure  level.  The  comparison  of  average  roof 
loadings  shows  that  the  predicted  curve  underestimates  the  initial  peak  load 
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and  overestimates  the  late  time  quasi-steady  loading.  Drag  coefficients  of 
-0.1  and  -0.3  are  more  correct  for  the  initial  load  and  the  late  time  loading, 
respectively.  A linear  relief  to  the  lower  pressure  level  in  a time  period 
of  1.5  transit  times  of  the  shock  across  the  structure  surface  would  adequately 
simulate  the  average  loading  behavior. 

The  experimental  and  predicted  overturning  moments  calculated  about  the 
lower  rear  edge  of  the  rectangle  are  presented  in  figure  38  for  both  overpres- 
sure levels.  These  moment  determinations  do  not  consider  the  resistance 
provided  by  the  rear  face  loadings  since  an  average  experimental  rear  face 
pressure  could  not  be  accurately  determined  with  only  two  measurements.  Never- 
theless, the  curves  may  be  used  with  the  average  front  face  and  roof  loadings 
to  provide  some  information  on  the  nature  and  degree  of  error  inherent  in  the 
simplified  techniques  for  determining  the  required  foundation  restraints  to 
prevent  overturning.  The  same  assumptions  concerning  rigid  body  behavior, 
inertial  effects,  and  structure  supports  as  made  for  the  arch  structure  are 
applied  to  the  rectangular  structure. 

Figure  39  shows  a comparison  of  the  experimental  and  predicted  vertical 
support  reactions  determined  for  the  100  and  300  psi  overpressure  levels.  A 
sign  convention  of  positive  values  representing  compressive  reactions  is  adopt- 
ed. At  both  overpressures , the  comparisons  indicatp  good  agreement  during  the 
diffraction  period  for  each  support  location.  The  later  time  comparisons  are 
not  quite  as  good  except  for  the  300  psi  downstream  support  where  excellent 
agreement  is  found.  The  100  psi  data  show  a relatively  large  overestimate  of 
compressive  reaction  and  underestimate  of  tensile  reaction.  An  underestimate 
of  tensile  reaction  is  also  observed  for  the  300  psi  data.  Both  the  over- 
estimate of  compressive  reaction  and  underestimates  of  tensile  reactions  are 
directly  related  to  the  late  time  overestimate  of  roof  loading  which  resulted 
from  the  prediction  technique. 

c.  Simulated  Ground  Plane 

When  surface  flush  elements  are  located  in  the  vicinity  of  abovegiound 
obstructions,  these  elements  may  be  subject  to  reflected  pressures  and  drag 
pressures  originating  at  the  obstruction.  Whether  these  effects  are  realized 
depends  on  the  location  of  the  element  with  respect  to  the  aboveground  protru- 
sion. The  "Design  Manual"  provides  a technique  for  estimating  the  loading 
on  surface  flush  elements  located  near  rectangular  aboveground  structures. 

This  technique  will  additionally  be  used  for  the  arch  structure  although  not 
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Figure  38.  Experimental  and  Predicted  Overturning  Moments 
Rectangular  Structure 
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Figure  39.  Experimental  and  Predicted  Vertical  Support  Reacti 
Rectangular  Structure  (Continued) 
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recommended.  The  "Design  Manual"  states  that  drag  effects  are  felt  within  a 
range  equal  to  the  height  of  the  structure  along  the  ground  surface  in  front 
of  or  behind  the  structure.  The  magnitude  of  this  drag  effect  is  taken  equal 
to  that  acting  on  the  front  or  rear  face  depending  on  where  the  surface  flush 
element  of  interest  is  located.  Reflected  pressures  are  assumed  to  act  over  a 
distance  KS  in  front  of  the  structure  where  S is  the  structure  height  and  K is 
a factor  depending  on  the  rarefaction  (CR)  and  reflected  (UR)  wave  speeds.  The 
value  of  K may  be  estimated  by  K = 2/ (CR/UR  - 1).  For  a distance  KS/2  in  front 
of  the  structure,  the  reflected  pressure  level  is  assumed  equal  to  that  acting 
on  the  front  face.  Over  the  remaining  KS/2  distance  the  peak  reflected  pres- 
sure above  the  incident  overpressure  level  is  assumed  to  attenuate  linearly  to 
zero  from  the  pressure  level  acting  on  the  front  face.  Allowing  D to  represent 
the  distance  of  the  point  of  interest  in  front  of  the  structure,  the  reflected 
pressure  is  assumed  to  vary  as  the  front  face  reflected  pressure  for  a period 
of  time  equal  to  |"| — |.  For  times  greater  than  this  value  or  distances 
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D greater  than  KS/2  the  reflected  pressure  level  is  assumed  to  attenuate  lin- 


early to  the  incident  overpressure  level  or  incident  overpressure  plus  drag 
pressure  if  D is  less  than  the  structure  height.  The  total 
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The  results  of  the  above  prediction  techniques  are  compared  with  the 
measured  data  in  figures  40  and  41  for  each  structural  shape  and  each  over- 
pressure level.  The  front  face  pressure  waveforms  (D=0)  used  for  the  predic- 
tions are  those  determined  earlier  with  the  techniques  presented  in  subsections 
4. a and  4.b. 


The  comparison  for  the  100  psi  arch  shown  in  figure  40a  immediately 
shows  that  the  prediction  technique  grossly  underestimates  the  drag  component 
of  loading  throughout  the  distance  equal  tL  the  structure  height.  The  compari- 
son additionally  shows  that  the  prediction  does  not  properly  account  for  the 
reflected  wave  behavior.  Near  the  structure  face  the  technique  results  in  an 
underestimate  of  reflected  pressure  level  because  the  predicted  average 
in!  ial  peak  pressure  level  across  the  entire  windward  face  was  used  in  the 
determination.  This  average  valjue  incorrectly  represents  the  reflection  level 
at  low  angles  of  incidence  near  the  ground  plane.  Away  from  the  structure  face 
{D  = 0.5S)  the  comparison  of  reflected  pressure  levels  is  more  favorable.  How- 
ever, the  improved  comparison  is  felt  to  be  only  coincidental.  At  the  most 
removed  point  (D  =.  1.0S)  the  data  suggest  that  the  reflected  wave  has  been 
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completely  eroded  by  rarefactions  and  that  drag  loading  effects  only  are  felt 
at  this  point  which  is  a reverse  of  the  prediction  technique.  The  later  time 
dynamic  pressure  effects  from  the  structure  are  seen  to  attenuate  with  distance 
from  the  structure  face.  This  attenuation  rate  is  much  lower  than  that  for  the 
reflected  wave.  The  indicated  attenuation  rate  is  such  that  the  dynamic  pres- 
sure effects  are  felt  at  a range  in  excess  of  the  structure  height.  The  under- 
estimation of  the  drag  loading  again  is  related  to  the  use  of  an  average  drag 
loading  across  the  windward  face  which  is  not  representative  of  the  pressure 
levels  near  the  ground  plane.  Had  the  predicted  reflected  pressure  been  de- 
termined using  a reflection  factor  for  zero  angle  of  incidence, the  near  struc- 
ture locations  would  have  shown  favorable  comparisons.  However,  the  locations 
farther  from  the  structure  would  have  shown  large  overestimates.  A similar 
relationship  would  develop  if  a zero  angle  of  incidence  drag  coefficient  was 

used. 

The  comparison  for  the  300  psi  arch  shown  in  figure  40  b again  shows 
an  underestimate  of  reflected  pressure  level  near  the  structure  face.  Unlike 
the  previous  comparison,  though,  the  data  show  very  good  agreement  between  pre- 
diction and  experiment  for*  locations  away  from  the  structure  face.  This  agree- 
ment also  extends  to  the  drag  loading  effects  occurring  in  later  time  at  the 
0.5  and  1.0S  locations.  An  immediate  observation  though  is  that  the  predicted 
levels  of  reflected  pressure  and  drag  loading  are  almost  equal.  This  observa- 
tion may  suggest  then  that  the  actual  quantity  being  measured  at  these  loca- 
tions is  drag  loading  only  and  that  the  reflected  wave  has  been  relieved  com- 
pletely by  rarefactions  within  a distance  less  than  one  half  the  structure 

height. 

The  experimental  data  for  the  arch  structure  have  shown  two  important 
factors  which  are  contrary  to  the  prediction  technique.  The  first  is  that  the 
reflected  pressure  level  attenuates  rapidly  with  propagation  distance.  The 
second  is  that  the  drag  loading  effects  are  not  constant  over  a distance  equal 
to  the  structure  height.  Although  not  proven  in  this  investigation  the  data 
additionally  suggest  that  drag  loading  effects  extend  past  a range  equal  to 
the  structure  height.  Because  of  the  limited  data  and  complexity  of  the 
loading  behavior  no  recommendations  for  changes  to  the  prediction  technique 
will  be  made.  Instead,  it  is  recommended  that  a preliminary,  analysis  use 
either  the  incident  overpressure  only  or  the  loading  at  the  lower  extremity 
of  the  structure  face  as  the  ground  surface  loading.  The  choice  of  which  to 
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use  will  be  dependent  on  the  type  of  structure  being  investigated  and  the  de- 
gree of  conservatism  desired.  Ideally,  additional  analysis  to  include  para- 
metric investigations  is  needed  which  will  define  the  attenuation  rates  of 
the  reflected  wave  and  drag  pressure  loadings  as  functions  of  structure  radius 
and  overpressure  level.  Once  this  is  accomplished  modification  of  the  pre- 
diction technique  is  warranted. 

The  comparisons  of  experimental  and  predicted  data  for  the  rectangular 
structure  are  presented  in  figure  41  for  each. of  the  overpressure  levels. 

These  comparisons  demonstrate  essentially  the  same  limitations  in  the  predic- 

■ 

tion  technique  as  were  found  for  the  arch  structure.  Specifically,  the  tech- 
nique  does  not  accurately  represent  the  attenuation  rates  of  the  reflected  wave 
and  drag  loading  effects.  Recognizing  the  similar  limitations,  the  same  recom- 
mendations concerning  preliminary  analyses  and  future  improvement  of  the  pre- 
diction technique  are  applicable  to  the  rectangular  structure.  Additional 
efforts  will  be  required  to  address  the  reason  for  the  double  reflection 
phenomena. 


SECTION  VI 
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SUMMARY 

The  objective  of  this  investigation  was  to  collect  data  describing  high 
pressure  shock  diffraction  associated  with  simple  arch  and  rectangular  struc- 
tural shapes.  Additional  information  was  also  collected  from  a stagnation 
probe  mounted  in  the  center  of  the  shock  tube  test  section.  The  data  collected 

were  then  used  to  determine  the  adequacy  of  existing  simplified  prediction 
techniques.  The  prediction  techniques  evaluated  were  those  proposed  by  the 
Air  Force  Manual  for  the  Design  and  Analysis  of  Hardened  Structures  (AFWL-TR- 
74-102).  In  satisfying  the  above  objective,  the  data  were  first  presented  in 
forms  which  provided  both  a general  description  of  the  loading  phenomena  and 
also  more  specific  characteristics  such  as  reflection  factors,  drag  coeffi- 
cients, and  clearing  times.  The  values  determined  for  these  more  specific 
characteristics  were  then  compared  to  values  proposed  by  the  preceding  "Design 
Manual."  Finally,  the  total  structural  loads  predicted  and  measured  were  com- 
pared so  that  the  accuracy  of  the  prediction  methods  could  be  observed.  The 
total  structural  loads  were  limited  to  average  structure  face  loadings,  over- 
turning moments,  structural  support  reactions,  and  ground  surface  loadings  in 
front  of  the  simple  structural  shapes. 

The  comparisons  between  recommended  and  experimentally  measured  reflection 
factors  resulted  in  the  following  findings.  For  a disc-shaped  plane  (stagnation 
probe)  the  recommended  values  of  reflection  factors  for  incident  overpressures 
between  90  and  330  psi  are  reliable.  The  arch  structure  tests  revealed  that 
reflection  factors  determined  on  an  inclined  planar  surface  at  various  angles 
of  incidence  may  be  used  to  represent  reflection  factors  along  the  arch  face 
at  similar  angles  of  incidence.  Naturally,  this  finding  is  applicable  only  to 
the  pressure  levels  investigated  (100  and  300  psi  overpressure).  The  major 
difference  between  the  recommended  and  measured  values  occurred  at  the  location 
of  the  extreme  angle  for  Mach  reflection  which  shifts  with  overpressure  level. 
The  recommended  curve  for  a planar  surface  indicates  a fixed  position.  The 
reflection  factors  determined  for  the  rectangular  structure  suggested  that  the 
recommended  values  are  somewhat  conservative  yielding  higher  values  at  the  100 
and  300  psi  overpressure  levels. 
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The  drag  coefficient  data  determined  during  the  stagnation  probe  tests 
revealed  a large  scatter  in  data  which  did  not  seem  to  be  a result  of  experi- 
mental error.  This  statement  is  made  based  on  the  repeatability  of  other 
measurements.  This  wide  data  scatter  suggests  that  dynamic  pressure  is  very 
sensitive  to  test  conditions  when  compared  to  the  other  measurements.  The 
comparison  between  the  median  variation  of  measured  drag  coefficient  and  rec- 
ommended value  as  functions  of  overpressure  level  showed  that  the  measured 
curve  indicated  a larger  increase  in  drag  coefficient  with  overpressure  levels. 
For  the  arch  structure  no  comparison  between  drag  coefficient  and  angle  of 
incidence  was  possible  since  no  recommended  curves  were  presented  for  the 
overpressures  of  interest. 

Comparisons  were  made,  however,  of  net  horizontal  drag  coefficient  using 
recommended  free  stream  cylinder  values.  The  comparison  showed  the  reconmended 
values  to  be  reliable.  Average  drag  coefficients  for  the  front  and  roof  sur- 
faces of  the  rectangular  structure  were  compared.  The  front  face  comparison 
showed  excellent  agreement.  Recommended  drag  coefficients  for  the  roof  struc- 
ture were  presented  only  for  overpressures  less  than  130  psi  which  allowed  only 
one  comparison.  The  recommended  value  for  the  100  psi  overpressure  test 
resulted  in  a large  overestimate  of  roof  surface  loading.  The  data  suggest 
that  an  average  drag  coefficient  for  the  roof  surface  should  be  a function  also 
of  the  structure  dimensions  and  not  just  overpressure  level. 


The  comparisons  of  total  structural  loading  for  the  arch  structure  showed 
that  the  prediction  techniques  significantly  underestimate  the  horizontal 
loading  during  the  diffraction  phase.  At  later  times,  the  techniques  are 
fairly  reliable.  For  the  average  vertical  loading,  excellent  agreement  was 
found  at  the  100  psi  level.  However,  a large  underestimate  of  vertical  loading 
was  found  at  300  psi  overpressure  after  the  diffraction  phase.  The  comparisons 
of  predicted  and  measured  overturning  moments  for  the  arch  structure  showed 
excellent  agreement  at  the  100  psi  level  and  poor  agreement  at  the  higher  level. 
This  poor  agreement  plus  the  erroneous  average  vertical  loading  in  the  post 
diffraction  period  resulted  in  underestimates  of  as  much  as  30  percent  for  the 
vertical  reactions.  A review  of  the  prediction  technique  to  determine  the 
reasons  for  these  differences  suggested  that  the  prediction  techniques  do  not 
adequately  address  reflection  phenomena  or  later  time  drag  pressure  variations 
with  structure  position.  Therefore,  instead  o*  recommending  changes  in  the 
input  values  for  the  techniques  it  would  be  more  appropriate  to  completely 
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revise  the  techniques.  Since  deficiencies  were  found  in  the  techniques  for 
gross  structural  behavior  it  is  not  recommended  that  they  be  used  for  refined 
analysis.  The  prediction  techniques  for  the  rectangular  structure  resulted 
in  excellent  loading  determinations  for  the  front  face  with  the  exception 
that  the  techniques  underestimated  the  clearing  time.  A revised  clearing 
time  parameter  is  recommended  which  includes  the  estimated  duration  of  che 
rarefaction  wave.  For  the  roof  surface  the  techniques  resulted  in  favorable 
loading  during  the  diffraction  process,  but  later  time  loads  were  significantly 
overestimated.  A post  diffraction  drag  coefficient  and  clearing  time  are  rec- 
ommended as  changes  to  the  prediction  techniques.  As  a direct  result  of  the 
overestimates  of  late  time  roof  loadings,  the  prediction  techniques  under- 
estimate the  overturning  moments  acting  on  the  structure  at  each  overpressure 
which  then  result  in  erroneous  vertical  structural  support  reactions.  The  com- 
parisons of  predicted  and  measured  ground  surface  loadings  showed  that  the 
prediction  techniques  are  seriously  in  error.  Revision  is  necessary  to  include 
changes  in  the  recommended  attenuation  rates  of  the  reflected  wave  and  attenua- 
tion rates  should  also  be  developed  for  the  upstream  effects  of  drag  loading. 
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APPENDIX  A 


MODEL  ST4  GAGE  SPECIFICATIONS 


Manufacturer 

Sensor 

Range 

Sensitivity 

Transducer  Capacitance 
Overload 

Maximum  Temperature 
Linearity 
Natural  Frequency 
Resistance 
Dimensions 


Susquehanna  Instruments 

Tourmaline 

10  - 10,000  psi 

0.1  pcmbs/psi 

10  pi co  Farads 

Maximum  Pressure  15,000  psi 

150°F  Gauge,  Intermittent  Gas  5,000° F 

± 2 percent  Full  Scale 

1.5  megahertz 

1010  ohms 

1/2  inch  diameter  by  2 inches  length 
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Table  B.l 

STAGNATION  PROBE  SHOTS 


INCIDENT 

MEAS.  MODEL  P™ 

NO. TOP  BOTTOM 

295  Stag. 

296  Probe  95  95 

297 

298 


Stag. 

Probe 


Stag.  135 
Probe 


Stag.  165 
Probe 


Stag.  167 
Probe 


Stag.  150 
Probe 


319 

320  Stag.  230 

321  Probe  : 

322 

323 

324  Stag.  320 

325  Probe 

326 


GAGE 

POSITION 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 

Incident 

Incident 

Stag. 

Stag. 


MIXTURES 
12  0„ 


100  Hi 


AMBIENT 

AMBIENT 

PRESSURE 

TEMP. 

(psi) 

(°c) 

12.01 

26.7 

12.02 


12.00 


11.99 


11.98 


12.09 


12.04 


12.04 


12.04 


AFWL-TR- 75-290 


Table  B.l 

STAGNATION  PROBE  SHOTS  (Continued) 


MEAS. 

NO. 

MODEL 

INCIDENT 

PRESSURE 

GAGE 

POSITION 

MIXTURES 

AMBIENT 

PRESSURE 

(PSI) 

AMBIENT 

TEMP. 

cc) 

TOP 

BOTTOM 

331 

332 

333 

334 

100 

100 

Incident 

Incident 

Stag. 

Stag. 

13  0, 
26  H‘ 

12.08 

26.1 

335 

336 

337 

338 

105 

103 

Incident 

Incident 

Stag. 

Stag. 

13  0, 
26 

12.06 

26.1 

339 

340 

341 

342 

■ 

110 

104 

Incident 

Incident 

Stag. 

Stag. 

13  0, 
26 

12.04 

26.1 

BTrS 

108 

100 

Incident 

Incident 

Stag. 

Stag. 

12.5  09 
25  H2  £ 

12.00 

24.9  ! 

347 

348 

349 

350 

105 

106 

Incident 

Incident 

Stag. 

Stag. 

12.5  09 
25  H2  6 

12.01 

25.0 

351 

352 

353 
.354 

E 

145* 

no 

Incident 

Incident 

Stag. 

Stag. 

12.5  0- 
25  H2  c 

12.00 

24.4 

355 

356 

357 

358 

ET® 

120 

105 

Incident 

Incident 

Stag. 

Stag. 

12.5  09 
25  H2  6 

11.97  j 

24.4 

359 

360 

361 

362 

no 

98 

Incident 

Incident 

Stag. 

Stag. 

12.5  09 
25  H2  6 

11  .98 

25.0 
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Table  B.2 

ARCH  MODEL  SHOTS 

MODEL 


INCIDENT 

PRESSURE 


BOTTOM 


105  100 


340  340 


105  100 


325  340 


110  100 


325  335 


105  100 


GAGE 

POSITION 


Incident 

Incident 

190° 

170° 

180° 

180° 


Incident 

Incident 

190° 

170° 

180° 

180° 


Incident 

Incident 

145° 

125° 

135° 

135° 


Incident 

Incident 

145° 

125° 

135° 


Incident 

Incident 

100° 

80° 

90° 

90° 


Incident 

Incident 

100° 

80° 

90° 


Incident 

Incident 

80° 

60° 

70° 

70° 


MIXTURES 


12.5  0, 


12.5  0, 
25  H.  L 


12.5  0„ 


12.5  0, 
25  Ho 


AMBIENT  I AMBIENT 
PRESSURE 


11.97  23.5 


11.97  23.0 


11.46  23.3 


11.94  23.6 


11.46  ?:.o 


11.97  21.5 


11.95  22.5 
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Table  B.2 

ARCH  MODEL  SHOTS  (Continued) 


INCIDENT 

PRESSURE 


GAGE 

POSITION 


Incident 

Incident 

10° 

350° 

0° 

0° 


Incident 

Incident 

30° 

10° 

20° 

20° 


Incident 

Incident 

50° 

30° 

40° 


Incident 
Inc! dent 
70° 
50° 


Incident 

Incident 

90° 

70° 

80° 


MIXTURES 


12.5  0, 


12.5  0, 


12.5  Or 


12.5  0, 


12.5  0, 


12.5  0, 


AMBIENT  I AMBIENT 
PRESSURE 
(PSI) 


11.94  22.8 


12.05  23.3 


12.06 


12.06 


12.04 


12.03  23.7 


12.02 


* Not  Recorded 


Mm 

u.  *;  . , • , . ; 


-Am. 
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Table  B.2 

ARCH  MODEL  SHOTS  (Continued) 

INCIDENT 


MEAS. 

PRESSURE 

GAGE 

NO. 

MODEL 

TOP 

BOTTOM 

POSITION 

81 

Arch 

Incident 

82 

Arch 

Incident 

83 

Arch 

105 

110 

170° 

84 

Arch 

150° 

85 

Arch 

160° 

86 

Arch 

160° 

87 

Arch 

Incident 

88 

Arch 

Incident 

8. 

Arch 

105 

no 

20° 

90 

Arch 

0° 

91 

Arch 

10° 

92 

Arch 

10° 

93 

Arch 

Incident 

95 

Arch 

310 

* 

0° 

96 

Arch 

350° 

97 

Arch 

10° 

98 

Arch 

10° 

99 

Arch 

Incident 

100 

Arch 

Incident 

101 

Arch 

285 

270 

30° 

102 

Arch 

20° 

103 

Arch 

40° 

104 

Arch 

Incident 

105 

Arch 

Incident 

106 

Arch 

300 

315 

50° 

107 

Arch 

40° 

108 

. Arch 

60° 

109 

Arch 

60° 

110 

Arch 

Incident 

111 

Arch 

Incident 

112 

Arch 

300 

310 

80° 

113 

Arch 

70° 

114 

Arch 

90° 

115 

Arch 

90° 

116 

Arch 

Incident 

117 

Arch 

Incident 

118 

Arch 

300 

235 

110° 

119 

Arch 

100° 

120 

Arch 

120° 

121 

Arch 

120° 

MIXTURES 


12.5  0, 


12.5  0, 


AMBIENT  AMBIENT 
PRESSURE  TEMP 

(psi)  (°c: 


12.01  23.7 


11.97  22.7 


12.00  23.7 


11.98  23.8 


11.96  24.0 


11.96  23.5 


11.96  24.5 


* Not  Recorded 


(touv 


fc  4 


i 

li 


290  275 


295  290 


305  295 


INCIDENT 

PRESSURE  GAGE 

TTOM  POSITION 


105  125 


310  325 


105  no 


295  290 


Incident 
Incident 
140* 
130° , 


Incident 

Incident 

170° 

160° 

180° 

180° 


Incident 

Incident 

20° 

10° 

30° 

30° 


70° 

80° 

90° 

90° 

Incident 

Incident 


70° 

80° 

90° 

90° 

Incident 

Incident 


70° 

80° 

90° 

90° 

Incident 

Incident 


70° 

80° 

90° 

90° 

Incident 

Incident 


MIXTURES 


12.5  0, 


12.5  0, 


AMBIENT  I AMBIENT 
PRESSURE 
(PSI) 


11.95  24.2 


11.95  24.3 


11.94  24.1 


12.04  20.5 


12.02  20.75 


11.92 


11.92  21.5 


* Not  Recorded 


fr. 


IMP*** 
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Table  B.2 

ARCH  MODEL  SHOTS  (Ccntinued) 


MEAS. 

NO. 


MODEL 


INCIDENT 

PRESSURE 


TOP]  BOTTOM 


GAGE 

POSITION 


Arch 

with 

wedge 

ex- 

ten- 

sion 


I 1.0  R 
0.5  R 
Stag. 
Stag. 
Incident 
Incident 


1.0  R 
0.5  R 
S 
S 

Incident 

Incident 


Arch 
wi  th 

wedge  280 
ex- 
ten- 
sion 


1.0  R 
0.5  R 
Stag. 
Stag. 
Incident 
Incident 


MIXTURES 

ambient 

pressure 

(PSI) 

AMBIENT 

TEMP. 

(°C) 

12.5  0o 
25  H2  2 

11.92 

20.5 

> 

12.5  0o 
25  H2  2 

11.91 

21.5 

44  09 
88  H2 

11.89 

22.0 

j 




mL.4'. 
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Table  B.3 

RECTANGULAR  MODEL  SHOTS 


MEAS.  MODEL 


PRESSURE 
topI  BOTTOM 


182 

183 

184 

185  Rec-  100 

186  tangle 


190 

191 

192 

193  Rec-  100 

194  tangle 


197 

198 

199  Rec-  315 

200  tangl e 

201 

202 

203 

204 

205 

206  Rec-  300 

207  tangle 


211 

212 

213  Rec-  100 

214  tangle 


217 

218 

219 

220  Rec-  280 

221  tangle 


* Not  Recorded 


GAGE 

POSITION 

MIXTURES 

AMBIENT 

PRESSURE 

(PSI) 

05F 
09  F 
11 F 
07F 
. 03F 
03  F 

Incident 

Incident 

12.5  09 
25  H2  £ 

11.97 

05F 
09F 
1 1 F 
07F 
03F 

Incident 

Incident 

12.5  09 
25  H2  c 

11.92 

05F 
09F 
1 1 F 
07  F 

Incident 

Incident 

44  09 
88  H2 

11.92 

05F 
09F 
1 IF 
07  F 
03F 
03F 

Incident 

Incident 

44  0, 
88  H2 

11.90 

05T 
09T 
1 IT 
03T 
03T 

Incident 

12.5  09 
25  H2  £ 

11.88 

05T 
09T 
1 IT 
03T 
03T 

Incident 

Incident 

44  09 
88  H2 

11.88 

AMBIENT 

TEMP. 

(°C) 


8 
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Table  B.3 

RECTANGULAR  MODEL  SHOTS  (Continued) 


. 


113 


"^11-11—1111  in  mu  i 
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Table  B.3 

RECTANGULAR  MODEL  SHOTS  (Continued) 


INCIDENT 

MODEL  — P-RE.SSURE — 
TOP  BOTTOM 

Rec-  110  170 

tangle 


GAGE 

POSITION 

05B 

Incident 

03B 

Incident 


05B 

Incident 

03B 

03B 


MIXTURES 


AMBIENT  AMBIENT 
PRESSURE  TEMP. 
(PSI)  (°C) 


12.5  0, 
25  H2  6 

12.03 

■ ■ 

19.75 

44  0, 
88  H2 

12.02 

20.5 

12.5  0, 
25  H2  1 

11.94 

21.1 

12.5  09 
25  H2  6 

11.93 

21.7 

12.5  09 
25  H2  *■ 

11.92 

23.0 

if  * 
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APPENDIX  C 


TEST  DATA  TRACES 


Table  C.l.  LOCATION  OF  COMPOSITE  DATA  TRACES 


Type  of  Data 


100  psi  Arch 


300  psi  Arch 


100  psi  Rectangle 


300  psi  Rectangle 


Inclusive  Page  Numbers 


116  through  179 


130  through  233 


234  through  272 


273  through  312 


?■ 

\*  l _ 
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4 ^ : * •» 


".vmrk  .T./IT wiyj 
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9|liH 
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mssmmm 
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MN  43,  44,  5C 


TIME(U-6CC 1 


T]M£(U-SECI 


AFWL-TR- 75-290 


800.-, 
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MN  49, 


lOO-Sfl-30 

56, 


i 

M 

I 

1 

| 

1 


PRESSUREl  PS  1 1 PR  £S3UR  1 P5 1 1 


PRESSURE l PS £ 


T1KECU-SEC) 


TIHE(U-SEC) 
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PRESSURE! PS [ I 


TIHE(U-SEC) 


TIHECU— SEC ) 
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T1MEIU-SEC) 
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PRESSURE!  PS  £ I PRESSURE! PS £ ) 
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PRESSURE ( PS 1 1 PRES SURE ( PS  £ I 


TIHE(U-SEC) 


lo.  100.  ISO.  200 

. 2 

SO.  300.  3! 
TIHE(U-SEC) 
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200.-1 

100-Sfl-INC*053 

«» vs*. 
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TINE(U-SEC) 


10O-5R-INOO6O 


TIMEOJ-SEC) 
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* 
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TIHE(U-SEC) 


100-SA-IN0157 


TIHE(U-SEC) 


PRESS*  REl  PSD  PRESSURE l PS  1 1 
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;i  v*- 
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175. -I 
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presslkeipsci  pressure!  psi 


m < mmrmm  W’  yvim  * 


TIHE(U-SEC) 
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TIHE(U-SEC) 


''Aiuyirpw-^ 
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TIHE(U-SEC) 
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1 ▼ I'm  ’ 

■gllllBSfi 
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Emmi 

PRESSURE ( PSr 1 PRESSURE! PS  I ) 


:ssuREiPsn  pressureipsh 


COMPOSITE  300-Sfl- INC 


TIHE(U-SEC) 
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PRESSURE! PSI ) PRESSURE! PS t J 
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TIHECU-SEC1 
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50.  100.  ISO.  200.  250.  300.  3S0.  400.  4S0.  500 

TIME(U-SEC) 


Li, 


- , . 


• a 


, v vv,  |,  ,*>«  . 

~v — -i-i 


(U-SEC1 


* W 

■HllyiaiKfir 

' t ii  t jftlL 

PSjB 

PJtl 

wwiiatRiii! 

|KTf 

mmm 

PRESSURE! PS I ) PRESSURE! PS [ ) 


-TT-m,  ■ -ns  l-‘TT>  '-.J'JI". 


■ 
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3200.- 


COMPOSITE  300-SR-30 
MN  101,  138, 


2800-- 


2H00  ■- 


2000- - 


1600.- 


200.  250.  300.  350. 
TIMECU-SEC ) 


sra^oF-5"-™ 


ETiirpwi^n 


500.  550.  600 


. V • 


i * 


«pr* 


COMPOSITE 


TIHEOJ-SECl 


TIHEOJ-SECl 


COMPOSITE  300-SR-90 
MN  11H,  160, 


TIMEUJ-SEC1 


COMPOSITE  300-SR-90 
MN  3T,  11H,  160, 


PRESSURE! PS t I 


COMPOSITE  300-SR-90 
MN  34,  114, 


TIHEOI-SEC) 


COMPOSITE  300-SR-100 
MN  32,  119, 


1^0.  2d0.  2SQ.  300.  3S0.  400-  450.  500 
TIHE(U-SEC) 


PRES5UREIP5 1 1 PRESSUREIPSI ) 
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TIHE1U-SECJ 
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TIHE(U-SEC) 


PRESSURE ( PST I PRESSURE l PS I 


T1MECU-SEC) 
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••  ■■■  ' v ' '■  - 
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400  q 

' 3Q0-Sfl-INC*007 

PRESSURE! PS i I 


^ _ 


PRESSUREl PS [ ) 
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TIKEOJ-SEC) 
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TIHE(U-SEC) 


•e "*v 
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iw. , •"»  * 


i*  -«*-  mffm 


TIME(U-SEC) 
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>0.  3067  ? 

TIHE(U-SEC) 


300-Sfl-  30*101 


TIHE(U-SEC) 
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300-Sfl-  20» 1 36 

TIMEOJ-SEC) 


T I ME  C U— SEC ) 
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TIHECU-SEC1 
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TIHE(U-SEC) 
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TIHE(U-SEC) 


PRESSURE ( PS t) 
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TIHE(U-SEC) 
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300-Sfl-  90*160 
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■'  s 
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TIHE(U-SEC) 
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" * * . 
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300-SA-190*009 
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miKKK 

TIHE(U-SEC) 
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PRESSURE (PS [I 


COMPOSITE  100-SR-05F 
MN  182,  190, 
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!0.  300.  3S0.  400.  4S0.  500.  550 
TINE(U-SEC) 
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100-SR-03T*214 


TIME(U-SEC) 


AFWL-TR-75-290 

000- -I 
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PRESSURE!  P51  i PRE5Sl*ElPSn 


TIHEtU-SEC) 


100-SR-11T*213 


TIHE(U-SEC) 
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3S0--1 

100-SR-Q9T*212 

TIHE(U-SEC) 


TIMECU— SEC ) 


AFWL-TR-75-290 
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100-SR-13T*234 

TIHECU-SEC1 
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100-5R-19T*232 
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TIMEfll-SEC) 
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400. -| 
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100-ER-15T*241 
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co 

a_ 

2 •" 
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U 
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Table  C.2 

LOCATION  OF  INDIVIDUAL  DATA  TRACES. 

STAGNATION  PROBE 

Gage 

Measurement 

Page 

Number 

Location 

Number 

Incident 

295 

316 

Incident 

296 

316 

FF 

297 

317 

FF 

298 

317 

Incident 

299 
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Incident 

300 
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FF 

301 
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FF 

302 
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Incident 
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FF 
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FF 
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FF 
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FF 
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FF 
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Tab  la  C..2  (Cont'd) 


Gage 
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Location 

Number 

Incident 
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